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Abstract 



These lecture notes provide a relatively self-contained introduction to field 
theoretic methods employed in the study of classical and quantum phase 
transitions. 



Chapter 1 
Introduction 



After the work of Wilson and others [T] in the 1970s on the renormahzation 
group (RG) in the theory of phase transitions, field theory methods originally 
employed in high-energy physics became an indispensable tool in theoretical 
condensed matter physics. Indeed, since the seminal work of Wilson, many 
excellent textbooks on field theoretic methods in condensed matter physics 
appeared, where a variety of other topics are also discussed; see, for exam- 
ple, Refs. [SlIallllElEllTllHlEllinillllin]. Onthe other side, the use of 
field theoretic methods in condensed matter physics led to further insights 
in particle physics. Thus, the study of classical phase transitions in several 
lattice spin models gave birth to lattice gauge theory (for a review, see Ref. 
[13]), which provided an approach to numerically tackle the strong-coupling 
regime of gauge theories. This interchange of ideas between these two fields 
became even more intensive after the discovery of high-temperature (high- 
Tc) superconductors, where it became apparent that traditional many-body 
techniques are insufficient to understand the mechanism behind high-Tc su- 
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per conductivity. For instance, gauge theories, both Abehan and non-Abehan 
were employed to understand several properties of doped Mott insulators 
[U [TJ [H], which eventually may become superconducting upon doping. 

These notes have the aim of introducing some field theoretic methods 
widely used in the study of classical and quantum phase transitions. Classi- 
cal phase transitions occur at a regime where quantum fluctuations do not 
play an important role, usually at high enough temperatures. In this case 
time, as obtained from quantum-mechanical equations of motion, does not 
play a role. Furthermore, in many situations is possible to use a continuum 
model Hamiltonian, which provides the stage for employing the field theory 
formalism. Wilson's RG was originally introduced in this context [Ij. The 
Hamiltonians are in this case actually determined by a Landau-Ginzburg 
type of expansion [T3] where the free energy is expanded in powers of the 
order parameter and its derivatives. The result is similar to Lagrangians of 
quantum field theories in Euclidean space, i.e., a space-time where time is 
imaginary |5]. It is this fact that allows for the application of field theory in 
the study of phase transitions ^21 13 El [TT] . Quantum phases transitions [6] , 
on the other hand, occur at zero temperature, such that time becomes impor- 
tant. The quantum phase transition is usually driven by some dimensionless 
parameter, like for example the ratio between two couplings appearing in a 
quantum Hamiltonian. Very often the Lagrangian for a system undergoing 
a quantum phase transition resembles one for a classical phase transition, 
excepts that one of directions of space is the (imaginary) time. Thus, when 
this occurs, the same field theoretic methods employed in the framework of 
classical phase transitions can be applied in the quantum phase. Sometimes 
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we say in such a context that a quantum system in ci + 1 spacetime dimen- 
sions (i.e., d spatial dimensions and one time dimension) is equivalent to 
a classical system having d + 1 spatial dimensions. Thus, if we think of a 
Landau-Ginzburg theory, we see that the effective Lagrangian of the quantum 
system has to look relativistic. There are several non-relativistic condensed 
matter systems possessing an effective description which is relativistic-like 
[31 El [71 Uni ttH [I2] • Here "effective" means a field theory description valid 
at some energy scale where the details of the underlying lattice model are 
irrelevant. There are cases, however, where time arises in the effective theory 
in a way which is not rotational invariant in d + 1 spacetime dimensions In 
those situations, frequency scales differently from momenta near the critical 
point, so they differ in scaling by some power P, |9]. Thus, we can write the 
scaling relation w ~ |p|^ at the critical point, which defines the so called dy- 
namic critical exponent, z. In this way, we see that a relativistic-like theory 
will have z = 1, so that the results from classical phase transitions can be 
applied to the quantum case more directly. We will see in these notes mostly 
examples of theories having z = 1. 

In these notes priority is given to the introduction of calculational meth- 
ods. Thus, the reader will usually find here very detailed calculations, which 
are often done step- by-step. 



When we say "rotational invariant" here, we are thinking of d+l EucUdean dimensions. 
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Chapter 2 



Ferromagnetism 

2.1 Spin in an external magnetic field 

Let us start with a very simple example spin dynamics, namely, the interac- 
tion of a spin with an external magnetic field. The Hamiltonian is 

if = -7S-B. (2.1) 
The Heisenberg equation of motion for the i-th spin component yields 

1^ = [^^'^] 

= -l[S.,S,]B,. (2.2) 

From the quantum mechanical commutation relation, 

^From now on we assume that repeated indices are summed over. 
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(2.3) 



we obtain immediately, 



dt 



(2.4) 



or, in vector notation. 



9S 



7(S X B). 



(2.5) 



The above equation simply describes the precession of a spin in an external 
magnetic field. 

The above equation can also be obtained from a Lagrangian for a classical 
spin. The classical spin is given by S = 5'n, where = 1. In this case we 
have the Lagrangian density. 



where the vector functional v4(n) has to be determined. This functional plays 
the role of a "momentum" canonically conjugate to n, so that ^(n) • 9(11 
would play the role analogous to a term pdq/dt in classical mechanics. In 
this context, we note that the Hamiltonian density T-L = —S'yn ■ B does not 
depend on A, which may appear to be a strange feature when compared to 
the most common situations in classical mechanics. However, we should not 
forget that wc have a constraint, = 1, which will play an important role 
in the derivation of the equations of motion. 



S[A{n)-dtn + -fn-B], 



(2.6) 
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The equation of motion for each component of the unit vector field n is 
given by the Euler-Lagrange equation: 



*(9(9tni) dui 
Thus, we have 9£/[9(9jnj)] = ^i(n), so that, 



(2.7) 



d{dtni) drij ^ 



We also have. 



orij V orii 



Altogether we obtain. 



dAi dA, 



(2. 



(2.9) 



(2.10) 



The term between parentheses is an antisymmetric second-rank tensor which 
is a function of the unit vector n. Therefore, it should have the form, 



dAi dAj 
drij drii 



(2.11) 



or, equivalently. 



OAk 



rii 



(2.12) 



The LHS of Eq. (2.12) is the i-th component of the curl of A. We can 



interpret ^4 as a vector potential defined in spin space. It is actually called 
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the Berry vector potential, and the corresponding term in the Lagrangian 
containing A, when associated to a functional integral representation of spin 
systems [6], constitutes the so called Berry phase. We will study the Berry 
vector potential further in a while. For the moment, let us proceed by deriv- 
ing the equation of motion for n. 



Using Eq. (2.12) in Eq. (2.10, we obtain easily, 



(2.13) 



Contracting both sides with eiminm yields. 



=^ {SijSmk - SimSjk)nmnkdtnj = 'yeimiUmBi 

=^ dm = 7(n X B),, (2.14) 

where from the second line to the third we have used both = 1 and 
n • dtn = (l/2)(9fn^ = 0. Therefore, we have just obtained the desired result, 
namely, 

dtn = 7(n x B) (2.15) 

Thus, once more a precessing vector is obtained, except that this time it is 
not an operator that is precessing. 
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2.2 The Landau-Lifshitz equation 



Eq. (2.15) is the simplest example of the so called Landau-Lifshitz (LL) 
equation. In order to study the dynamics of ferromagnetic systems, we have 
to go beyond the situation of an external magnetic field. To do this, first 
note that 



15H , , 



where 



H = -jS j ciVn- B. (2.17) 

If H has a more complicate functional dependence on n, involving spatial 
variations of the unit vector, B will be itself a function of n, constituting 
in this way an effective magnetic field. By assuming spatial isotropy, such a 
Hamiltonian is given by 



n= (fr 



JS 



2 



-din ■ diU — --fSn ■ B 



(2.18) 



2 

where JS"^, with J > 0, is the bare spin stiffness. The first term in the 
Hamiltonian above can be motivated via the continuum limit of the lattice 
spin Hamiltonian for a ferromagnet, 



H = -J^Si-Sj, (2.19) 

(id) 
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where the lattice spin fields are assumed to be classical vectors, Sj = Srij, 
and the lattice sites are summed over nearest neighbors. Indeed, by inserting 
the lattice Fourier transformation 



V L 



(2.20) 



where L is the number of lattice sites, in the Hamiltonian above, we obtain. 



H = Y,J{^) Sk-S_k, 



(2.21) 



where 



(2.22) 



a=l 



and we have set the lattice space equal to unity. In the continuum limit, we 
have, 



j^(k)^-6J+il 



(2.23) 



The second term above contributes to the first term in (2.18), while the term 



—6 J just adds an irrelevant constant to the Hamiltonian, since 



- 6 J ^ Sk • S_k = -6 JS^ 5Z = -GJS^L. 



(2.24) 



Let us consider now the LL equation with an effective field determined 



by the Hamiltonian (2.18). The LL becomes 
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dtn = SJ{n X V^n) + 7(11 x B). (2.25) 

Next we set B = and study small transverse fluctuations around the 
e3-axis, i.e., 



n=^e3 + 5nx, (2.26) 
V2 

where 



Sn± = ^[cos(k • r — ujt)ei + sin(k • r — ujt)e2]. (2.27) 
v2 

This solution describes a so called spin-wave, whose quanta are known as 
magnons. By inserting the above spin-wave profile in the LL equation, we 
see that it solves it provided 



J 02 

u = -^e, (2.28) 
which yields the spin- wave spectrum for a ferromagnet. 



2.3 The Berry vector potential, magnetic monopoles 
in spin space, and hedgehogs in real space 

Let us discuss the Berry vector potential further. If we compute the flux of 
the curl (in spin space) of ^(n) through the unit sphere 5^2 (the bar reminds 



us that 5*2 is embedded in spin space) using Eq. (2.12), we obtain 
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dSi e 



S2 



dSi Hi 

S2 

2" ,^ r , fdn dn 



(2.29) 



where we have used the parametrization n = (sin 77 cos sin sin ^, cos 77), 
and e Z is the winding number. This result shows that the Berry vector 
potential corresponds to the field strength of a magnetic monopole in spin 
space. 

Note that this result is supposed to be an intrinsic property of the Berry 
vector potential such that should hold beyond the case of a single spin in 
an external field. However, for the case where many spins are involved and 
have a gradient energy, we have a further topological aspect. In fact, when 
n has a nonzero gradient, we have the following fiux through a sphere S2 
in real space (note the absence of the bar over 5*2 and the differences in the 
integration measure). 



'S2 



r dxi dxm 
dip J deeum-Q^-^eijun ■ {djU x ^^11) 

I dip J d9{5ij5rr,k - ^i'^^j^^~Q;^~^^ ' ^^i^ ^ ^''^^ 
dxi dn dx„i dn 

X 



(2.30) 



Jo ^'^ Jo V dxi ' ~ de dx, 

/ dip j d9n 
Jo Jo 



/ dn dn 
\d^ ^ d9 
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We see that the last hne of (2.30) is twice the second hne of (2.29). Thus, 
we have, 



dSieijkn ■ {djU x ^n) = 47riV. (2.31) 

Therefore, we see that there is not only a topological object living in spin 
space, but there is also one living in real (three-dimensional) space, a so called 
hedgehog. Just like the magnetic monopole in spin space, the hedgehog in real 
space is also a point-like topological object. These topological objects arise 
also in two-dimensional quantum antiferromagnets, which will be discussed 
later in these notes. 




2.4 The classical nonlinear a model 

In order to study the static critical behavior of a ferromagnet, we can consider 



just the Hamiltonian (2.18) and try to compute the partition function via 



the functional integral representation, 

Z = J Vn6{n^ - 1)6"^ (2.32) 

where the S function enforces the constraint. We can rescale the coordinates 
and the temperature such as to have only 1/T as the coefficient of the gradient 
term H and remove the prefactor JS"^. We will also generalize the model such 
that it will have d dimensions and the unit vector field will have n components 
instead of three. In this way, we will be able to consider two approximation 
schemes, namely, one where an expansion in e = c? — 2 is made, and another 
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one where the number of components n is large, such that an expansion in 
1/n can be done. 

2.4.1 The expansion in e = d — 2 

The main idea behind all e-expansions is to make a perturbative analysis 
of the problem around the dimensionality where the coupling constant of 
the model is dimensionless. In such at the so called critical dimen- 

sion (i.e., the dimension that makes the coupling constant dimensionless), 
perturbation theory is relatively well behaved. The perturbation series is 
still divergent and frequently needs to be resummed, especially when higher 
orders in perturbation theory is involved. However, rigorous mathematical 
results exist in some cases where out of some (initially) perturbation series a 
non-perturbative expansion may be rigorously constructed [16]. The situa- 
tion is considerably more difficult if we are not at the critical dimension and 
the theory approaches the critical point, thus becoming strongly coupled. In 
these cases there are several resummation procedures; for a review on these, 
with focus on variational perturbation theory, see the textbook by Kleinert 
and Schulte-Frohlinde [8]. 

For the classical non-linear a model, the coupling constant is the tem- 
perature T. Simple dimensional analysis shows that this coupling becomes 
dimensionless at d = 2. We will show later that for d = 2 the model is 
actually asymptotically free. We are ultimately interested at = 3, so the 
goal will be to set e = 1 at the end of the calculations. The case d = 2 and 
n = 2 is special and will be treated separately in the next chapter. 



13 



The calculations will be performed up to second order in the temperature 
and will serve as a good introduction to the calculation of loop integrals 
occurring in Feynman diagrams. Therefore, we will perform all integrals 
exactly in d dimensions and in great detail below. Simple properties of 
analytical continuation of the dimension in the integrals will be assumed. 
This is known as dimensional regularization. The properties of dimensional 
regularization are very simple and are discussed in several textbooks; see for 
example Ref. |H] 

By resolving the constraint = o"^ + tt^ = 1, the Hamiltonian of the 
nonlinear a model can be written as 



H 



1 

2T 
1 

2r 



1 -7r2 



(2.33) 



where tt = (tti, . . . ,7r„_i). 

The Green function can be written as [n = (cr, tt)] 



G{x) = (n(x)-n(O)) 

= ((T(x)a(O)) + (7r(x) ■ 7r(0)) 

= (Vl - 7r2(x) - 7r2(0)) + {Tzix) ■ 77(0)). (2.34) 
If one rescales tt by ^/T and expand up to order T, we obtain 
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G{x) = l-|(7r2(x) + 77^(0)) +T(7r(x)- 77(0)) + 0(T2) 

= l + (n-l)T[G'o(x)-Go(0)] + 0(T2), (2.35) 



where 



Go(x) = / — , (2.36) 



and we have used the short-hand notation 



(2.37) 



(27r)'^' 

Note that we have used translation invariance to write {-K^ix)) = {^^(0)) = 
(n-l)Go(O). 

The dimensional regularization rules demand that (see, for example, the 
textbook [8] and Appendix |B|) 



^ 0, for all a, (2.38) 



a 



q 



m 

which immediately leads to 6*0(0) = 0. 

In order to calculate Gq{x) explicitly we use Feynman parameters: 

Go{x) = dX e^P-^-^P' (2.39) 
Jo Jp 

By performing the Gaussian integral in p (by just completing the squares), 
we obtain 
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Goix) = ^^°°rfAA-^/^exp (-g) . (2.40) 
The substitution u = a;^/(4A) yields 



r)d—2 poo 

2'='-2r(d/2 - 1) 1 



(47r)d/2|a.|d-2 {d-2)Sd\x\'^-^'- 



(2.41) 



where Sd = 2tx'^/'^ /V{d/2) is the surface of the unit sphere in d dimensions 
(see Appendix |A]) and we have made use of the definition of the gamma 
function: 

T{z) = / dTT'-^e~\ (2.42) 
Jo 

By performing the substitution r = as, with a constant, the above definition 
of the gamma function gives another very useful formula in calculations with 
dimensional regularization: 



1 



dss'-^e-"'. (2.43) 







r{z) 

To illustrate the usefulness of this formula, let us calculate two functions that 
will be needed later: 



and 
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Up) 



g2[(p + g)2]-/2 



(2.45) 



where a G 



Using Eq. (2.43), we obtain for Eq. (2.44), 



T(a/2)yo n 



(4vr)'^/2r( 



4r y 



(4vr)d/2r(a/2)|x|'^-° 



2-1 -u 



(4vr)'^/2r(a/2)|x| 



d—a ' 



(2.46) 



Note that by setting a = 2 we recover the expression for Gq{x). 



The calculation of laip) is more involved. By using Eq. (2.43), we can 
write 



-1 pea POO f Jd„ 

'"^^^ - r(a/2) X Vo 'J W ^ ^ 



After performing the Gaussian integral in g, we obtain 



Up) 



oo /-oo 



{4n)^/^T{a/2) Jo Vo (n + r2)'^/' 



(iri / dT2 



a/2-1 



e.p,-^),(2.4S) 



Now we perform the change of variables 
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ri = ra, = (1 - r)a, (2.49) 

where < r < 1 and < a < oo. Note that the Jacobian of the transforma- 
tion is (J. The resuh is 



1 /■! /"OO g-p2r(l-r)o- 

After introducing the change of variables s = j9^r(l — r)cr, the above integral 
can be rewritten as 



\ \d-a-2 pi (d-a-2)/2 poo 



(47r)<^/2r(a/2)" 



r f ^±^1 r f ,,/(^),,(?.6i) 



(47r)'^/2r(a/2) V 2 / V 2 )T{d-l-a/2) 
where from the second line to third we have used the integral 



'fa"(l-.)'=r(l+^Ml±^. (2.52) 







r(2 + a + h) 

We can easily go up to third order in T using the methods illustrated in 
the paper of Amit and Kotliar [T7]. Up to second order in T, for example, 
we can proceed as before and expand further a/1 — Ttz'^[x) to obtain 
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(Vl-T7r2(a:)Vl-T7r2(0)) = I - {n - l)TGo{Q) - —{'k\x)tt\Q)) + . . . 



= l-^(7r2(a;)7r2(0)) + ..., (2.53) 



T2 



where we have used dimensional regularization to set Go(0) = 0. The remaing 
avarage involving four pi fields can be decoupled with the help of Wick's 
theorem, i.e., 



{Tz'ixWm = 2[(7r(x).7r(0))-(7r^(x))(7r^(0))] 

= 2{n-l)[Gl{x)-Glm='^{n-l)Gl{x). (2.54) 



The third order contribution is more difficult to obtain, since in this case 
the non-linear factor 



in the partition function has to be taken into account when performing the 
correlation function average. By expanding the exponential factor above in 
powers of T we obtain also derivative contributions. The latter are more 
easily handled with in momentum space. Afterwards we can transform the 
result to real space. Thus, we have the following expansion. 




(2.55) 
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{2txY5\p) + T(n - 1)^ + -T\n - l)h{p) 



T\n-l) 



where [T7] 



(2.56) 



and 



h{p) 



(2.57) 



/c(p) 



(p + qY 



JqJqi q^ql{p + q + qif' 



(2.58) 



J ( I f {p + q?{p + qi? 

'^^^ P'jjq.q'ql{p + q + qi?' 



(2.59) 



(2.60) 



J n J 



{q - qi) 



(2.61) 



-g^gi q^ip + qYqlip + qiY' 
Note that the integral h appearing in the term of second order in T is 



just the Fourier transform of Eq. (2.54), taking into account also Eq. (2.53). 



However, it is instructive as one more exercise in dimensional regularization 
to perform the integral h- This can be done straightforwardly by considering 
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Eq. (2.51) with a = 2 to obtain 



h{p)=cid)\p\ 



d-A 



(2.62) 



where 



c{d) 



T{2- d/2)T'^{d/2-l] 
{A-Kyi^Vid - 2) 



Using Eq. (2.46), we obtain 



(2.63) 



h{x) = Gl{x), 



(2.64) 



which corresponds to the expected result, in view of Eqs. (2.53) and (2.54) 



Integrating over qi in Jc, using the result of J;,, we obtain 



Up) 



c{d) 



By using Eq. (2.51) with a = —d, we obtain 



(2.65) 



Up) 



cid)T{l-d)Tid/2-l)Tid) (,_3) 
{A7rY/2T{-d/2)T{3d/2 - 1) 



Next we can use Eq. (2.46) to obtain Ic{x): 



(2.66) 



(2.67) 



The integral Id{p) can be written in the form 
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Id{p) = Jd{p)/p' + Kd{p)/p^ + p-^ 



{p + qf 



q^{p + q + qiY'' 



(2.68) 



{p + qf 



q J q\ 



q^qlip + q + qiY' 



(2.69) 



and 



q J qi 



2p ■ qi{p + qf 
q^ql{p + q + qiY' 



(2.70) 



The last integral in Eq. (2.68) vanishes due to the rule (2.38) (by integrating 
over qi). 



Using (2.51) with a = 2 to integrate over gi, we obtain 



Jd = c{d) 



g2[(p + g)2](2-d)/2 



(2.71) 



Now we use once more (2.51), but this time with a = 2 — ci, to obtain 



Up) ^ Jd{p)/p' 



r(2 - d/2)T%d/2 - l)T{d - 2)T{d - 1) 
{AnYT{d - 2)r(l - d/2)T{3d/2 - 2) 



|p|2(d-3)_ (2.72) 



For Kd we need first to calculate 



{P + g)^^ 



qif. 



Ql 



qfip + q + qiY 



(2.73) 



I.e., 
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{p + 
1 



gi • (p + g) 

(p + q + qiY - {p + qf 



ql 



2{p + qY,,^ 
-1 / I 



qlip + q + qiY 

c{d) 



[{p + 



2l(rf-4)/2 



(2.74) 



where rule (2.38) has been used in the first and last terms of the second line 



of the above equation. Using this result in the expression for Kd-, we obtain 



KM = -c{d) 



p-{p + q) 



g2[(p + ^)2]{2-d)/2 



-p^Jd - c{d)Ld{p), 



(2.75) 



where 



Ld{p) 



p-q 



,g2[(p + g)2](2-d)/2 

1 f{p + qf 



p2 _ q,2 



2 g2[(p + g)2](2-d)/2 
P'^Icjp) _ P^Jdjp) 

2c{d) 2c{d) 



(2.76) 



At the end, we obtain 



23 



Since hip) = V^Iliv)) we have 

Je(x) = C\d) 



|p|2(3-d)' 



which can be solved using (2.46) with a = 2(3 — d). We have 



iJx 



c-(rf)23(^-2)r(3d/2-3) 



(47r)'^/2r(3-d)|x|3(rf-2)' 
For If{p), it is easy to see by expanding {q — qi)^ that 



The first integral above vanishes because of (2.38), and we obtain 



q^^ 



We can now use (2.74) with p + q replaced by p to obtain 



.^^M|p|2(rf-3)_ 



Therefore, 
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c\d)2^('^-')T{3d/2 - 3) 
^ ~ 2(47r)'^/2r(3-d)|x|3('^-2)- ^^-^"^^ 

In the expression for G{x) it is the following combination that appears: 



(2.84) 



Therefore, 



= 1 + T(n- 1)6-0(0;) + y(n-l)G'^(x) 



(n - l)r='G3(x) 



(n + l)ci - 4 (n - 2) cos(7rd/2)r(2 - ci/2)r(3rf/2 - 3)" 



6(3rf-4) 2r(rf-2) 
+ C(T^). (2.85) 

The above is the final expression for G{x) up to third order in T. By 
setting d = 2 + e and expanding for small e, we obtain exactly the same 
expression as in the paper of Amit and Kotliar [1^, i.e., 



G{x) = l + T{n-l)Goix) + —{n-l)Gl{x) 



+ 



r3(n- 1) 

{n - 2y 
4 



n — 3 {n — 2)e 



6 



+ 4(^-2) 



6 



C(3) 



(2.86) 



Note that due to the expression for Go{x), Eq. (2.86) contains poles for e = 0. 
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These poles can be subtracted by introducing two renormalization constants 
flEl [T7]. Z and Zt, such that the renormahzed correlation function reads 



Grix;t) = Z-^G{x-T = Zt^i'H), 



(2.87) 



having no poles for e = 0. 

The bare correlation function should be independent of the renormaliza- 
tion scale /i, which is expressed by the equation, 



or, using the chain rule. 



(2. 



/^^ + /3W^ + CW G{x;T = Ztfi-H) = 0, (2.89) 
where the renormalization group (RG) functions P{t) and ({t) are given by 



/3(t) = fi 



dt 



at) = /i 



d\nZ 



(2.90) 



Using the derived perturbative expansion for the correlation function, we 
obtain. 



/3(t) =et-{n- 2)t^ - (n - 2)t^ 



(2.91) 



at) = {n-l)t + -{n-l){n~2)t\ 



(2.92) 



The vanishing of the /3 function determines the critical temperature as an 
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expansion in powers of e, 



n-2 



i- 









C(e^). (2.93) 



In terms of the RG functions above, the exponents u and rj are respectively 
given by [17] 



1 



V = Cite) - e, 



corresponding to the critical behavior 



(2.94) 



where ^ is the correlation length, and 

G{x) ~ 



— V 
C) 1 



1 



\X 



d-2+r] ' 



(2.95) 



(2.96) 



The correspondence between the RG functions and the critical exponents 
follows by making dimensionless variables in G{x) explicit. For instance, we 
can write. 



where A is the ultraviolet cutoff. Thus, we have the behavior. 

Therefore, near the critical point we have, 
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(2.97) 



(2.98) 



dlnZ 
9/i 



(2.99) 



Comparison with Eq. (2.92) yields the second of Eqs. (2.94) 



For the critical exponent u, we linearize the /3 function near the critical 
point tc, i.e., 



d{t - 1, 



-/3'(g(t-tc), 



(2.100) 



where the prime denote a derivative with respect to t and the minus sign 
reflects the negative slope of the /3 function at tc (ultraviolet stability of the 
fixed point). Integrating the above equation yields 



A 



(t - Q 



(2.101) 



which upon the identification /i = leads to the determination of the 
critical exponent v. 
Explicitly, we have 



1 



e + 



n-2 



+ 0(a, 



(2.102) 



n-2 



1 - 



n — 1 
n-2 



e + 



n{n — 1) 

2(72-2)2' 



(2.103) 



Note that the calculations allow to obtain r] with one order of e higher than 
for 1/z/. 
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2.4.2 The 1/n expansion 
The saddle-point approximation 

A common nonperturbative approach to study the non-hnear a is the large 
n limit [6l |5] , where n is taken to be large while ng is kept fixed. Such an 
analysis is more easily done if the constraint is implemented with the help of 
a Lagrange multiplier field, such that the Lagrangian becomes 



£=^[(9,n)2 + ^A(n2-l)]. (2.104) 



This Lagrange multiplier field arises from the partition function (2.32) when 



the functional integral representation of the delta function is used: 

5(n2 - 1) = y I5Ae-^^'^'"'^("'-^). (2.105) 

Let us integrate n — 1 components of n exactly and call the non- integrated 
one cr. The resulting effective action reads 



S,^ = i!i_ilTrln(-V2 + i\) + ^ j rf'^x[a(-V' + i\)a - ^A]. (2.106) 

The limit n — )■ cxd is obtained from the saddle-point approximation to the 
above effective action. We will consider a uniform saddle-point with i\ = 
and a = s. Thus, from 



= 0, ^ = 0, (2.107) 

om'^ as 
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we obtain, 



= 1 -nT / ^— , (2.108) 



m^s = 0. (2.109) 
For T <TcWe have s 7^ 0, so that = 0, leading to 

s' = l-nT [ (2.110) 
J (27r)<ip2 ^ ^ 

At the critical point, T = Tc, we have that s = 0, so that 



d'^p 1 SdA 



d-2 



nT^ J (27r)'^p2 (27r)<^(ci - 2) ' 
where A is the ultraviolet cutoff. Therefore, 



(2.111) 



s'^^iT,-T), (2.112) 

^ c 



which implies that the critical exponent of the order parameter is /3 = 1/2. 
For T > Tc, on the other hand, we have s = and 7^ 0, such that 



r d-^p 1 _ 1 
J (27r)'^p2 + ^2 - (2.113) 



If we now write 



1= t( f ^ f d'^p 1 f d'^p 1\ 

\J {2'KYp' + m? J (27r)<ip2+y (277)'^^^^' 
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we obtain, 



T-Tr 



nTm 



(2.115) 



The integral above can be evaluated with help of the integral Ji of Appendix 
|B| Indeed, it is proportional to the integral Ji evaluated in c? — 2 dimensions. 



I.e., 



{2'kYp^{p^ + m? 



1 ^ , , , 2^-'^'K-'^l^m'^-^ ^ f d 

— -hid - 2 = r 2 - - 

27rc;-2 ^ ' d-2 \ 2 



Thus, 



(2.116) 



— r 2 - - = -. 

d-2 \ 2 T' 



or, after substituting the value of Tc, Eq. (2.111) 



(2.117) 



m = A 



d/2-2 



T{d/2)T{2-d/2) 



1 - 



{27ry{d-2)A^-'^ 
SduT 



l/{d-2) 



(2.118) 



It turns out that the mass gap m = ^ ^, where ^ is the correlation length. 
Therefore, the corresponding critical exponent is 



d-2 

The limit d ^ 2 of Eq. ( |2.118[ ) yields 



(2.119) 



^We are making use of the well-known result = lim„j_s.oo(l + x/mY 



31 



A / 27r\ , , 

m=-exp -— , 2.120 
7T \ ni J 

which imphes that for d = 2 the system is gapped for all T > 0, i.e., no 
phase transition occurs in this case. This is consistent with the Mermin- 
Wagner theorem [20], which states that a continuous symmetry associated 
to quadratic dispersions cannot be broken in d = 2. Note that the Mermin- 
Wagner theorem rules out the symmetry breaking at ci = 2, but not neces- 
sarily the phase transition. Indeed, we will see in the next chapter that d = 2 
and n = 2 is special. In fact, although no symmetry breaking occurs in this 
case, a phase transition happens even in its absence. 

The P function for the dimensionless coupling t = A'^^^T at large n is 
easily obtained by demanding the scale invariance of the mass gap, i.e.. 



. dm 

A- = 0, (2.121) 



which yields, 

and we see that for d = 2 the theory is asymptotically free. 
1/n corrections to the saddle-point 

Now we want to compute the 1/n corrections to the saddle-point calculation. 
Let us consider for instance the theory at the critical temperature and find 
the anomalous dimension of the a-field. This is achieved by first expanding 
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the effective action (2.106) up to quadratic order in A, such as to obtain the 
A-propagator. This yields 



5". 



n 



eff 



d'^x J d'^x'Goix-x')Goix'-x)\{x)\{x') + ^ J d'^xa{-V^+i\)a, 

(2.123) 



where Go{x) is given in Eq. (2.36). Thus, it is easy to obtain the correlation 



function for the field A in momentum space as 



(A(p)A(-p)) 



{2/n) 



f Jl2 
J (2n 



d'ig 1 



(2.124) 



or, using the result (2.62), 



(A(p)A(-p)) 



c{d)n\p\ 



d-i ■ 



(2.125) 



where c{d) is given in Eq. (2.63). In order to find the anomalous dimension 



of the (T-field we need to compute the cr-propagator up to order 1/n. This 



is done by taking into account the vertex iAcr^/(2T) in Eq. (2.123). The 
lowest order contribution to the a-propagator features two vertices. The 



corresponding Feynman diagram is shown in Fig. 2.1 Thus 



cr 



{p)a{-p)) 



-1 



p^ + 



p^ + 



d'^q {X{p + q)X{-p - q)) 
(2^ ^2 

d'^q 1 



nc{d) 



{2ttY \p + q\'^-^q'^' 



(2.126) 
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Figure 2.1: Feynman diagram representing the 1/n correction to a- 
propagator. The internal sohd hne corresponds to the free cr-propagator, 
while the dashed line is the A-propagator. 



The above result contains an integral which can be evaluated with the help of 



the integral (2.45) along with its explicit evaluation in Eq. (2.51). However 



for a = — 4, like in Eq. (2.126), Eq. (2.51) has a pole, which is actually 



related to a logarithmic behavior. The trick to complete the calculation is 



to replace d — 4 in Eq. (2.126) by a, evaluate the integral explicitly via Eq. 



(2.51 ), and perform an expansion in powers of d — 4 — a, setting a = d — A 



at the end. 



The dangerous contribution in Eq. (2.51) comes from the factor r((2 + 
a — (i)/2), whose singularity for a = — 4 can be isolated as 



2 + a-d 



+ (regular terms). 



(2.127) 



2 J d-A-a 

Furthermore, the result of the integration is proportional to jp]'^"""^, which 
can be expanded as 



^d-4-ap2g(d-4-a)ln(p/A) ^ X'^-i-»pi[l + (rf _ 4 _ c^) ln(p/A) + ...]. (2.128) 

Now we note that the anomalous dimension is given by the infrared behavior 
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of the propagator, i.e., 

HpM-p)) ~ (2.129) 

Thus, before setting a — d — 4, we neglect the terms which are smalUer in 
comparison with p^~^ as p — > 0, anticipating aheady that r) > 0. Thus, by 
keeping only the dominant terms in the infrared, we can safely set a = o? — 4 
to obtain 



{a{p)a{-p)) ^ ^ 
which should be compared to 



2(d-4)r(d-2) /p- 
nr(2-d/2)r(d/2 + l)r2(d/2-l) VA. 



For d — Z this yields 



(2.130) 



p^-'' -77ln(p/A)], (2.131) 
to obtain, after some simplifications. 



_ {A-d){d-2fnd-2) 
^~ ndT{2-d/2)T^{d/2) ' ^ ^ ' 



V-ir^- (2.133) 

The above result for 77 is the same as obtained for the 0{ri) Landau-Ginzburg 
model at large n. Indeed, both models have in the framework of a 1/n 
expansion the same critical exponents, belonging henceforth to the same 
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universality class [5] 



In order to see the non-perturbative character of Eq. (2.132), let us 



compare this expression with the perturbative expansion in powers of e 



d — 2 in Eq. (2.103) for n large. If we expand Eq. (2.132) in powers of e, we 
obtain 



n \ 2 



(2.134) 



which is precisely Eq. (2.103) in the limit n ^ 1. The large n result agrees 



with the large n regime of the e-expansion even up to order four. Indeed, 
we have that up to this order the anomalous dimension is given by (see for 
example Ref. [5j and references therein) 



rj = g+(ri-l)e2|-l + -g 



+ 



^ I 2 — n 3 — ^ , , 
-6+(n-2) + C(3) 



i'\ + 0{e'), (2.135) 



where we have defined 



n-2' 



(2.136) 



and 



h = -^{n' - 22n + 34) + ^C(3)(n - 3). 



(2.137) 



In the limit 1 Eq. (2.135) becomes 
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,2 

n \~ ' ' 2 



= - ( 1 - e + ^ - i±Me3 ) + 0(6% (2.138) 



which agrees precisely with the expansion of Eq. (2.132) in powers of e up 
to the same order. This comparison clearly exhibits the non-perturbative 
character of the 1/n expansion, since the lowest non-trivial 1/n correction to 
1] is already able to reproduce the large n limit of perturbation theory up to 
fourth order. 
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Chapter 3 

The Kosterhtz-Thouless phase 
transition 

3.1 The XY model 

When n = 2, the local constraint n| = 1 in the classical Heisenberg model is 
solved by writing 

rij = (cos6'i,sin6'i), (3.1) 
such that the Heisenberg Hamiltonian becomes 

H^-JS^Yl ^^^i^i - ^j)- (3-2) 
{hi) 

Such a planar magnetic system is known as XY model. As we will see in a 
later Chapter, this model is in the same universality class as a superfluid. 
The reason is not difficult to understand, since universality classes are usually 
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determined by the underlying symmetry group, which in the present case 
is 0(2), and this is the same as the unitary group U{1). The spontaneous 
symmetry breaking of the U{1) group will be studied in several ways for d > 2 
in the next Chapters. Here we will explore the case n = 2 and d = 2, where 
no spontaneous symmetry breaking occurs, as we have already mentioned in 
the last Chapter. While for > 2 no phase transition can occur at = 2, 
the situation is different for n = 2. When n = 2 and ci = 2 a phase transition 
occurs without spontaneous symmetry breaking. 

3.2 Spin- wave theory 

Let us consider the classical non-linear a model for n = 2. By writing 
n = (cos ^, sin 0), we obtain. 



-H = ^(V^)^ (3.3) 



The above equation is just the continuum limit of the XY model in Eq. (3.2 ), 
up to a trivial renaming of the couplings. 

It is easy to see that the correlation function G{x) = (n(x) • n(0)) of the 
non-linear a model can be written in this case as 

G{x) = (e^[^(-)-^(o)]). (3.4) 

More explicitly. 



""^Note that translation invariance implies (cos 6'(a;) sin 0(0)) = (sin 6'(x) cos 0(0)) . 
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G{x) = pee-^'^'^1^(^'^)'-''(^')'^(^')], (3.5) 
where Z is the partition function and 

J{x')=i[6\x-x')-6\x')\. (3.6) 

Note that we are not setting d = 2 yet. We will see soon why it is convenient 
to do so. We can perform the Gaussian integral over 9 exactly to obtain 



G{x) 



(fx' / dV'J(x')e?(x'-x")J(x") 



exp 

exp{r[^(x)-^(0)]} 



(3.7) 



where 



g{x) 



(3. 



We are not going to use dimensional regularization in this Chapter. It is 
actually essential to us not to have ^(0) = 0. This is important, in order 
to take the limit (i — )■ 2 in a clean way. Thus, by evaluating ^(0) explicitly 
using a cutoff, we obtain. 



(27r)^(d-2)' 



(3.9) 



From Eq. (2.41), we obtain 
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Thus, 



which has a straightforward hmit c/ — )■ 2, 



C2 



-(A|x|) 



2-d 



(3.11) 



g{x)-g{o) = ~^HA\x\). (3.12) 



Therefore, we have 



(A|x|) 

where 



= (3.13) 



= (3.14) 

is the anomalous dimension of the theory. Interestingly, in contrast with 
the case d > 2, the correlation function for n = 2 and d = 2 features a 
temperature-dependent anomalous dimension. This is a particularity of the 
two-dimensional XY model, as shown first by Kosterlitz and Thouless [2T] . 

The above result is exact within spin-wave theory, which corresponds to 
a situation where the fact that 6 is an angle is not taken into account. The 
periodicity of 6 is, however, crucial for characterizing the phase structure of 
the model. This will be the subject of the next Section. 
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Before closing this Section, let us point out the agreement between the 
perturbation theory for the non-linear a model for n = 2 and the results of 
this Section. Indeed, if we set n = 2 in Eq. (2.85), we obtain. 



G{x) = 1 + TG,{x) + ^G,{x) + ^Gl{x) + O(T^). (3.15) 

2 

The above expansion contains the first terms of the expansion of exp[TG'o(3;)]- 
Note that Gq{x) is the same as G{x) in this Section. Furthermore, since we 
are not using dimensional regularization, we must make the replacement 

Go{x)^g{x)-g{o). 

Another point worth mentioning in the context of the perturbation theory 



of the previous chapter is that for n = 2 and d = 2 the (3 function (2.91) 
vanishes, i.e., 

dt 

Thus, for n = 2 and d = 2 perturbation theory is just reproducing the result 
of spin- wave theory, i.e., a free theory. 



3.3 Two-dimensional vortices 

We have already mentioned in the previous Section that the spin-wave theory 
neglects the periodicity of 6. Thus, the spin-wave analysis implies 



dx-Ve = 0. (3.17) 



c 
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The above is not always true in the case of a periodic field pTl [2], |22] . This 
is precisely the case if we interpret as the superfiuid velocity In the 
presence of vortices, the circulation of the superfiuid velocity is quantized, 

^ dx-V9 = 27in, (3.18) 

where n G Z is the winding the number counting how many times the 
closed curve C goes around the vortex. The equation above is just the 
Bohr-Sommerfeld quantization condition applied to a superfiuid. 

In two dimensions vortices are just points. Thus, the problem of studying 
the statistical mechanics of many vortices in two dimensions is much easier 
as in three dimensions. In three dimensions vortices are either infinite lines 
or loops. The statitiscal mechanics of vortex loops is much more complicated 
and requires often the use of special duality techniques involving multivalued 
fields PES]. 

Let us give a concrete example of a single vortex at the origin in two 
dimensions. Such a vortex has to be singular at x = (xi,X2) = (0,0). Such 
a vortex configuration is simply given by 

V^= ^(-X2,xi). (3.19) 

Thus, 

dx.Ve = ^-^^^^^P^. (3.20) 
x^ 

^The superfiuid velocity is given in terms of the phase of the order parameter as [32] 
V, = {h/m)Ve. 
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Using polar coordinates, 



Xi = rcos^, X2 = 'rsin^, (3.21) 

where r = we obtain simply, 



/.27r 

dx-\/e= de = 27r, (3.22) 
c Jo 



which is just Eq. (3.18) with n = 1. It is easy to see that the 9 corresponding 



to such a singular configuration is given by 

e = arctan (^^-^ . (3.23) 

We are of course interested in a configuration involving many vortices. 
We will consider only vortices with vorticity n = ±1, which are energetically 
more favorable. A many-vortex configuration can be easily constructed. We 
have. 



VOv = 5^ 7^r^^[-(a; - ^^)2el + {x- Xi)ie2], (3.24) 



where Xi is the position of the i-th vortex and qi = ±1. Note that V6v looks 
like an electrostatic field caused by a potential By in two dimensions. In this 
analogy, the vorticities qi play the role of point charges in two dimensions. 
Thus, we are going to study the statistical mechanics of a two-dimensional 
Coulomb gas. Since 



Vln(A|a;-a;o|) = T^^, (3.25) 
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we have that the vortex contribution of the Hamiltonian is given by 



= ^Y^QiQi j d^xVHK\x-Xi\)-V\n{K\x-Xj\). (3.26) 
A partial integration yields 



Hv = ^Y.'li'l^ j d^^H^\x - Xi\) ■ [-V'ln(A|x - (3.27) 

and the boundary term vanishes by imposing the "neutrality" of the two- 
dimensional Coulomb gas, 



J]gi = 0. (3.28) 



Since 



VMn|x-Xo| =27r5^(x-Xo), (3.29) 



we have, 



Hv ^ -^^qiqjln{h\xi- Xj\). (3.30) 

The electric susceptibility of this Coulomb gas is obtained by coupling 
the vortex Hamiltonian to an external "electric" field and taking the second 
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derivative of the free energy with respect to it. The resuh for a single dipole 
is a susceptibihty proportional to (x^). In calculating this average the Boltz- 
mann factor involving the energy of a vortex-antivortex pair has to be used 
as weight for the averaging. Thus, 

{x^) ~ / rfrr^e-^-/^)!'^^^^). (3.31) 

Note that the integral features a short-distance cutoff. The integral converges 
only if 



T < ^. (3.32) 



We obtain in this case, 



(x^) -J . (3.33) 



Therefore, the susceptibility will diverge at the critical temperature 

Tc = |. (3.34) 

At this critical temperature the system changes the phase in a singular way 
(the susceptibility diverges). For T < Tc the system is in a dielectric phase 
of vortices. For T > Tc the vortices unbind and we have a plasma of vortex 
"charges" ±1. This phase transition between vortex states is the celebrated 
Kosterlitz-Thouless phase transition |21j . 

At the critical temperature the anomalous dimension has the universal 
value 
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viTc) = \. (3.35) 



3.4 The renormalization group for the Coulomb 

gas 

In this Section we will derive the RG equations governing the KT phase 
transition. These equations were derived for the first time by Kosterlitz [23] 
and are important in order to derive one of the most remarkable results of 
the KT phase transition, namely, the existence of a universal jump in the 
superfluid density at [21] • This result was later confirmed by experiments 

In this Section we will derive the RG equations for the Coulomb gas us- 
ing a scale-dependent Debye-Hiickel screening theory. The derivation will 
be done in d dimensions, as the three-dimensional result will be useful to us 
later on in a different context. The RG equations for a d-dimensional Coulomb 
gas were first derived by Kosterlitz [26j using the so called "poor man scal- 
ing" [27]. The Debye-Hiickel method used here to analyze the d-dimensional 
Coulomb gas follows Refs. [2HI 122] , which is inspired from a paper by Young 
[30] . who analyzed the two-dimensional case. 

To begin with, let us consider the bare Coulomb potential in d dimensions. 



which can be written using the result (3.11) as 



f/o(r) = -An'KoVir), (3.36) 
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where we have defined r = |a;| and Kq = 1/T, and 



V{r) = g{r)-g{o) 

{d-2){2nY 



2-d 



(3.37) 



with a = A ^. From the bare Coulomb interaction (3.36) we obtain the bare 
electric field, 



(3.38) 



dr SdT^ ^ 

The crucial step for our analysis is the introduction of a scale-dependent 
dielectric function defining an effective medium for the Coulomb system, 
so that the renormalized electric field determines a renormalized Coulomb 
potential f/(r), i.e., 



E(r) 



dU 
dr 



(3.39) 



dr Sd£{r)r'^^'^ 

In order to establish a selfconsistent equation, we need to specify the dielectric 
constant via the electric susceptibility of the system. Thus, from the standard 
theory of electricity, we know that 



e{r) = l + SdX{r). 



(3.40) 



where the scale-dependent electric susceptibility is given by 
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r 

d~l. 



X{r) = Sd dss''~'a{s)n{s), (3.41) 

J a 

and the polarizability for small separation of a dipole pair is 

a{r) ^ (3.42) 

(Ji 

while the density of dipoles is given by a Boltzmann distribution in terms of 
the renormalized Coulomb potential, 

n{r) = z^oe-^^'\ (3.43) 



with zq being the bare fugacity. Therefore, by integrating Eq. (3.39) the 



selfconsistent equation for the renormalized Coulomb potential is obtained: 



Now we set I = ln(r/a) and define the effective coupling via 



K-\l) = '±^e^^-'^^. (3.45) 

Since 



dr Saeiae^yi'^-^) ' ^ ' 

we obtain 
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Furthermore, 



dl Kq dl 

(Jj 

= {d-2)K-^ + z^, (3.48) 



where z'^{l) is obviously defined by the second term of the second fine in the 
equation above. Thus, the renormahzed fugacity z{l) satisfies the differential 
equation. 

By introducing the dimensionless couplings k = a^^'^K and y = a'^z, we 
finally obtain the desired RG equations for the d-dimensional Coulomb gas. 



^ = id-2)n-' + y\ (3.50) 
dl 



^d-^4^]y. (3.51) 



dl V 5*^ 

For d = 2 the above equations yield the celebrated RG equations for the KT 
phase transition 
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^ = {2-7rK)y. (3.53) 



Eq. (3.53) has a fixed point at Kc = 2/tt, which corresponds precisely to 



the critical temperature = 7r/2 obtained before. However, Eqs. (3.52) and 



(3.53) contain additional information. The fiow diagram actually features a 
line of fixed points for k > Kc at zero fugacity. Usually we say that the KT 
fiow diagram has a fixed line rather than a fixed point. The RG fiow diagram 



is shown in Fig. |3.1[ Note that the we set 2 — ttk as the horizontal axis, such 
that the critical point occurs for a vanishing abscisse. The blue line shown 
in the figure is a separatrix delimitating three different regimes. Below the 
separatrix and for 2 — vr/t < we have a dielectric phase of vortices, which 
corresponds to the low temperature regime where a vortex-antivortex pair is 
tightly bound, forming in this way a dipole. The dielectric phase is gapless 
because dipoles do not screen [31]. That is the reason why for 2 — ttk < 
there is a fixed line. At a fixed point all modes are gapless, so if there is a line 
of fixed points, like in the KT case, the excitation spectrum is completely 
gapless along this line. For 2 — ttk > 0, on the other hand, Debye screening 
occurs and a gap arises. This is a (classical) metallic phase (or plasma phase) 
of vortices. 

Let us study in more detail the RG equations. To this end we introduce 
the new variables, 



X = 2-7iK, Y = -^y. (3.54) 



The RG equations become 
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y 




2-JTK 
Figure 3.1: Schematic flow diagram for tlie KT pliase transition. 



f (3.55) 



f=.Yr. (3.56) 

Tlius, tlie family of liyperbola 

- f2 = const (3.57) 
are RG invariants. We introduce further new variables, 

u = X + Y, (3.58) 

and 

v = Y -X, (3.59) 
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such that the RG invariance equation becomes 



u{l)v{l) = UqVo 



(3.60) 



where Uq = u{0) and Vq = v{0) are given initial conditions. Thus, 



du 



-{uv + u^) 



(3.61) 



It is straightforward to solve the above equation by direct integration. The 
result is 



Now we set r = ^, i.e., we let the distance scale be equal to the correlation 
length. In this case we have that the mass gap is approximately given by 



which is a result characteristic of the KT transition. Note that we do not 
obtain in this case a power law for the mass gap. 

For d > 2 the situation is completely different. First of all, the Coulomb 
gas cannot be interpreted as vortices any longer, since in three dimensions 
vortices are one-dimensional objects, lines or loops [2]. However, there are 




(3.62) 




(3.63) 
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physical systems in three dimensions with point-hke topological defects where 
an analysis similar to the one made here is applicable. For example, there 
are systems where magnetic monopole-like defects occur in three spacetime 
dimensions [28l|29l|32l|33l[Sl[35l[Sll371|38llM^ Second, for > 2 



the RG equations (3.50) and (3.51 ) do not have nontrivial fixed points. Thus, 
no phase transition occurs in this case. The excitation spectrum is always 
gapped, so that the system remains permanently in the plasma phase. 
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Chapter 4 



Bose-Einstein condensation and 



superfluidity 



4.1 Bose-Einstein condensation in an ideal gas 

The Lagrangian for an ideal Bose gas is written in the imaginary time for- 
mahsm as 

C = b*drb- ii\b\^ + —Vb* -Vb. (4.1) 
2m 

All thermodynamic properties of the ideal Bose gas can be derived from the 
partition function, which is given by the functional integral representation, 



Z = I Vb*Vb e"^, (4.2) 

where 
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S = J dr J d'^rC. (4.3) 

The functional integral above is to be solved using periodic boundary condi- 
tions 6(0) = 6(/3) and 6*(0) = h*{l3). 

We see from the action for the ideal Bose gas that the case /i = is 
special. Indeed, if /i = 0, the action is invariant by a transformation where 
the Bose field is shifted by a constant, h ^ h -\- c. Note that the periodic 
boundary conditions make the contribution c*drh vanish. Thus, the special 
role of the /x = regime can be accounted for by shifting the Bose field by a 
constant, i.e., 

b = bo + b, (4.4) 

and we require that bo minimizes the action. This requirement implies that 
no term linear in b or b* appears in the action. This is only true provided 



libo = 0. (4.5) 

This equation is fulfilled either for bo — and 7^ 0, or 6o 7^ and /j, — 0. 
The Lagrangian is rewritten as 

C = -fi\bo\^ + b*(^dr-ix- ^V'^ b, (4.6) 

where the Laplacian term is obtained through partial integration in the ac- 
tion. By performing the Gaussian functional integral over b, we obtain, up 
to a constant, the result. 
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exi)(;)'l'//|6o|-) ^^^^ 



where V is the (infinite) volume. Thus, the free energy density is given by, 



/ = ^InZ 



|2 I ■'■ i„ j„j- fa .. ^ V72 



-M' + ^ Indet ( a. - /. - — ) . (4.8) 



Since the determinant of an operator is given by the product of the eigenval- 
ues of the operator, we have to solve the differential equation, 

5r-/^-^V2^^ = E^, (4.9) 

where E is the eigenvalue. The equation above should be solved with periodic 
boundary conditions ip{0) = In order to solve the eigenvalue problem 

we perform a Fourier transformation in the spatial variables. 



In this way the partial differential equation becomes an ordinary differential 
equation of first order. 



2 

dr-l^+^] i'ir, p) = E^(t, p), (4.11) 



2m 

which can be easily solved to obtain. 
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V'(r,p) = ^/'(0,p)exp 



T[E + fl 



p 

2m 



(4.12) 



Due to the periodic boundary condition, the above equation for r = /3 be- 



comes 



1 = exp 



/3 + 



2m 



(4.13) 



This imphes, 



^n(p) 



-iUn - /i + 



2m' 



(4.14) 



where a;„ = 2?™ / /3 with n G Z is the so called Matsubara frequency. Inserting 



these eigenvalues in Eq. (4.8) yields, 



/ = -H'-oP + ^ / 1. (-.. - A + l^j . (4.15) 

Our interest is to compute the particle density, n, which is the variable 
conjugated to the chemical potential. We have. 



n 



which yields. 



(4.16) 



df-p 



(2vr) 



(4.17) 



In order to have 6o 7^ we need /i = 0, so that the above equation becomes 
for 6o 7^ 0, 
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(2^ 



(4.18) 



lUJr. 



2m 



The Matsubara sum appearing above is performed in the Appendix D Thus 



n 



\bo\' + 



cfip 



1 



(4.19) 



e 2m — 1 

The remaining integral is over a Bose distribution for free bosons in d di- 
mensions. An integral involving a more general spectrum is evaluated in 



Appendix [Cj Using Eq. (C.9) with z = 2 and c = l/(2m) and making some 
simplifications, we obtain, 



n 



1 - 



ad/2) fmT\ 



n 



d/2 



(4.20) 



Note that we have solved for |&oP) which is the so called condensate density. 
Its physical meaning is that for fi = the particle density zero momentum 
gets depleted due to temperature effects. The density at zero momentum 
emerges because in momentum and frequency space. 



b{Un, p) = bo6'^{p)6nfl + b{Un, p) 



(4.21) 



such that Bq is associated to the zero momentum and zero Matsubara mode 
contribution of the Bose field. The Bose-Einstein condensation is thus the 
macroscopic occupation of the zero momentum state, in which case b repre- 
sents the fluctuation around the condensate. 



Eq. (4.20) can be rewritten as 
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= n 



1 




(4.22) 



where 



2ti { n 



1 



m [ad/2) 



(4.23) 



is the critical temperature. For T = the condensate vanishes. Note that 
for T = all the particles are condensed. This is a feature of the ideal Bose 
gas. We will see that in the interacting case the condensate is also depleted 
at T = due to the interaction. 

From the expression for the critical temperature we see that it vanishes 
for d = 2, implying that no condensate exists in a two-dimensional ideal Bose 
gas at finite temperature. This result is actually more general and holds even 
in the interacting case. It is known as Hohenberg's theorem |19j . 

4.2 The dilute Bose gas in the large limit 

In Chapter [2] we have studied the 0{n) classical non-linear a model in the 
large n limit. We will now use this knowledge to perform a expansion 
for an interacting Bose gas. Such an expansion actually corresponds to the so 
called random phase approximation (RPA) for the dilute Bose gas introduced 
long time ago [IHl EOl E21 E31 El] ■ That the expansion for the dilute Bose 
gas corresponds to RPA was recognized by Kondor and Szepfalusy [53] long 
time ago. Their analysis will be revisited here from a functional integral 
point of view |I3] . 
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4.2.1 The saddle-point approximation 

Let us consider the following action for a A'"-component interacting Bose gas: 



a=l ^ ' \a=l / 



, (4.24) 



where and 6* are complex commuting fields. The partition function is 
then given by 



-,-S 



(4.25) 



In order to perform the 1/A^-expansion we introduce an auxiliary field A(r, r) 
via a Hubbard-Stratonovich transformation: 



N , 
a=l ^ 



2m 



(4.26) 



Now we integrate out N — 1 Bose fields to obtain the effective action 



= {N-l)Tr\n{dr-pi- — + iX 



dr / d\ 



,.|a.-,-^ + a)6+lA^' 



, (4.27) 



where we have called b the unintegrated Bose field. 

Next we extremize the action according to the saddle-point approxima- 
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tion (SPA), a procedure that becomes exact for N ^ oo. This part of the 
calculation is practically identical with the one for the ideal Bose gas. This 
is done by making the replacement iX — )■ Aq and b — )■ bo, with Aq and 60 be- 
ing constant fields, followed by extremization with respect to these constant 
background fields. From this SPA we obtain the equations 



(Ao -/i)6o = 0, 



(4.28) 



Ao = 9\bo 



Ng 



1 



(2it) 



lUJr, 



(4.29) 



The large A^ limit is taken with Ng fixed. Below the critical temperature 



Tc we have 60 7^ 0, and thus from Eq. (4.28) Aq = yU. Therefore, Eq. (4.29) 
becomes 



provided d > 2. The second term on the RHS of thee above equation is, up to 



the prefactor A^, the same as the one in Eq. (4.20) for the condensate density 
of the ideal Bose gas. The particle density is obtained as usual n = —df/dfi, 
where f = —InZ/ {NVP) is the free energy density. This gives us 



n 



M 

A^ 



d'^p 
(27rp 



exp 



1 

/5 (1^ + ^-/^ 



(4.31) 



By setting Aq = in Eq. (4.31) and using Eq. (4.30), we obtain 
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n 



and therefore the condensate density becomes, 



no 



TV 



n 



d/2 



where 



(4.32) 



(4.33) 



2% 

m 



n 



2/d 



(4.34) 



We see that the SPA does not change the value of Tc with respect to the 
non-interacting Bose gas. Indeed, the SPA corresponds to the Hartree ap- 
proximation and it is well known that it gives a zero shift. 



4.2.2 Gaussian fluctuations around the saddle-point ap- 
proximation: Bogoliubov theory and beyond 

In order to integrate out A approximately, we consider the 1/A^-corrections 
to the SPA by computing the fluctuations around the constant background 
fields 6o and Aq. By setting 



h = hQ + h^ i\ = \q + iX, 



(4.35) 



and expanding the effective action (4.27) up to quadratic order in the A field, 
we obtain 



63 



h*[^r-^^ + Xo~j^]h + ilipih + h,h* + \h\^) + 



■y /^dr f dr' I d\ I dVA(r,r)G'o(r-r',r-r')G'o(r -r,r'-r)A(r',r'), (4.36) 



^0 



where S^^^ is the effective action (4.27) in the SPA and 



^o(^'r) = -pY. j (^^^^""""^"^^ol^cun,?), (4.37) 



with 



G'o(ia;n,p) 



iw„ + /i - Ao - - 



2 • 

2m 



(4.38) 



After integrating out A the effective action can be cast in the form 



^eff = Sll"^ + -Tr In [6{t - t')6\v - r') - Ng Go{r -r',r- r')Go(r' - r, r' - r)] 

^(r, r)r(r - r , r - r')^"^(r', r')^(r , r ) 



+ 





+ ^ j^^ dr j d'^r j^^ dr' j d'^r' ^^t^v) 







1 



r(r-r',r-r')^^(r',r')^(r',r') 



+ \j^dT j d'^r J^t' j rfV^t(^^r)^(r,r)r(r-r',r-r')^"^(r',r')^(r',r'04.39) 



1 




where we have introduced the two-component fields 
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¥{t, r) 



6*(r,r) 6(r,r) 



6(r,r) 



(4.40) 



which satisfy = 2|6p. The matrix M(r — r',r — r') has a Fourier 

transform given by 



M{iuJn,p) 



-iUn + S{Un, p) bl t{iUn, p) 

(6^)^f(iw„,p) + £:(u;„,p) 



where 



+ Xo + ^ + \bo\'^t{iUn, p) 
2m 



with 



(4.41) 



(4.42) 



t{iujn,p) 



9 



(4.43) 



1- NgU{lLOn,p) 

which is the Fourier transform of the effective interaction r(r — r',r — r' 
and 



i]-{iun, P) = ^ X] / 7^-^'^o(«c^n + «c^m, P + q)G'o(«a;m, q), (4.44) 



is the polarization bubble. The effective interaction can be represented in 



terms of Feynman diagrams as in Fig. 4.1 Physically A corresponds to the 
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+ 



+ 



+ 



Figure 4.1: Feynman diagram representation of the effective interaction Eq. 



(4.43). The dashed hne represents the bare A-field propagator while the 



double dashed line represents the dressed 1/A^-corrected A-field propagator. 
Cont inuou s lines represent 6-fields and each loop is the polarization bubble 
Eq. (4.44) formed by two 6-field propagators in convolution. The effective 



interaction is obtained as a geometric series of polarization bubbles. 



fluctuation of the particle density and thus the effective interaction (4.43) 



gives in fact the density-density correlation function. An effective interaction 



like the one in Eq. (4.43) was already obtained some time ago by a number 



of authors [IHl EOl EB [S2] • Thus, the 1/A^-expansion is actually equivalent to 
a random phase approximation (RPA) considered previously in the literature 



IHl |52]. Explicit evaluation of the Matsubara sum in Eq. (4.44) yields 



n(icu„,p) 



(p + q) 

2m 



2m 



+ Ao - /i 



(4.45) 



where nB{x) = l/(e^^ — 1) is the Bose distribution function. 
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By inverting the matrix (5.31) we obtain the propagator 



(4.46) 



where 



2 „ ' 

1^ + Ao - /i + |6oPr(ia;n, p) - |6o|^r2(ia;„, p) 



(4.47) 



and 



Ffe.fo. (iWn,P) 



(6*)2r(zu;„,p) 



2 

^ + Ao - /i + |6oPr(ia;„, p) - |6o|^r2(ia;„, p) 

(4.48) 



is the anomalous propagator. 

For T > TcWe have that &o = and the matrix propagator becomes diag- 
onaL In this case if we neglect the interaction terms from the effective action 



(4.39), the b propagator corresponds to the Hartree approximation. Thus, 



above Tr we have to consider the effective interaction between the bosons in 



Eq. (4.39) in order to obtain a nontrivial result for the excitation spectrum. 



This is achieved by computing the 1/A^-correction to the propagator. Below 
Tc, however, a nontrivial result for the excitation spectrum is obtained from 
the pole of the propagator even without considering the correction to 



it. This is easily seen from the structure of the propagators (4.47) and (4.48) 
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where the effective interaction f (zci;„,p) appears exphcitly. Above Tc the ef- 
fective interaction appears in the propagator only in the next to the leading 
order in 1/N. 

4.2.3 The excitation spectrum below Tc 

As we have discussed in the previous Subsection, below a nontrivial re- 
sult for the excitation spectrum is obtained in an approximation where the 
interaction term of the effective action is neglected. Thus, we now undertake 
a study of the spectrum of the system in such a Gaussian approximation. 
Later we shall see that the effective action in the Gaussian approximation 
gives the free energy density up to the order 



From the pole of the matrix propagator (4.46) we obtain that the energy 
spectrum E{p) satisfy the equation 



E\p) = Re||- + Ao-/i+|6o|'f[^(p)+^'5,p] 1^ - |6o|'Re f ^[^(p) + ^5, p] 

2 \ 2 / 2 

|^ + Ao-/ij +2(^ + \o-(x]\bo\'ReT[Eip)+t5,p], (4.49) 



where 6 — i- 0^. Note that Eq. (4.49) can be written as the product of two 



elementary excitations, -E^(p) = Ei{p)Et{p), where 



Ei{p) = ^ + \o-ix + 2\bo\^Re T[E{p) + i5, p] 



(4.50) 



Et{p) = ^ + Ao - /i, 
2m 



(4.51) 
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are the spectrum of the longitudinal and transverse modes, respectively. 
When Xo = fi we obtain that the transverse mode is gapless, consistent 
with Goldstone theorem. 

By inserting the saddle-point value Aq = /i, we obtain the following self- 
consistent equation for the excitation spectrum 



where the notation E{p) = £'(p)|ao=^ is used. Note that the above spectrum 
corresponds to a generalization of the well known Bogoliubov spectrum |18] . 
The difference lies in the fact that in the 1/A^-expansion the coupling con- 
stant g is replaced by the effective interaction T[E{p) + iS,p] |l9l |50]. At 
zero temperature Il{iu, p) vanishes and the excitation spectrum corresponds 
to the usual Bogoliubov spectrum. As we shall see, the modification of the 
spectrum by the effective interaction accounts for thermal fluctuation effects 
in higher temperatures and allows us a consistent treatment of critical fluc- 
tuations near Tc. 

In Eq. (4.52) we can legitimately replace |6oP by Nn, since the error 



committed in such a replacement is of higher order in 1/N. Thus, Eq. (4.52) 
becomes 



^(P) = \l^2+ — R«r[E(p) + p]. (4.53) 



Note that since T[E{p) + i5,p\ ~ C(l/A^), Eq. ( |4.53[ ) is independent of N 
for —7- oo. 

In order to obtain the spectrum of elementary excitations we need to 
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evaluate the polarization bubble (4.45). Unfortunately, it cannot be evalu- 



ated exactly, although many of its properties and asymptotic limits can be 
worked out exactly [521 [55]. For instance, for distances much larger than the 
thermal wavelength the polarization bubble can be evaluated exactly |52l |55] . 
This is called the classical limit in the early literature of the field. In the 



classical limit we can approximate the Bose distribution in Eq. (4.45) by 
nB{x) ~ l/(3x. In such a limit we can write 

n(iw„,p) = no(p) + ni(iw„,p), (4.54) 

where 



no(p) = -4m^T j 



{2t,Y [(p + q)2 + 2m(Ao - /i)] [q^ + 2m(Ao - /i)] ' 

(4.55) 



IIi(zC(;„,p) = QuvTiuj. 



(^.TiY 2miu)n - p2 - 2p ■ q 

^ (4.56) 



[(p + q)2 + 2m(Ao - /i)][q2 + 2m(Ao - /i)] ' 
Setting Aq = /i, we obtain the following result for the effective interaction: 



where p = |p| and 
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V{2-d/2)V\d/2-l) 

When 2 < (i < 4 we obtain for small p that fo(p) ~ p'^~'^ / {adTm?N), such 
that the excitation spectrum is given approximately by 



E{p) - J-^ P^'~'^'\ (4.59) 

The above equation reflects the dynamic scaling behavior, E{j)) ~ p^ of the 
excitation spectrum, where z is the dynamic critical exponent. Thus, from 



Eq. (4.59) we see that it implies a dynamic exponent 



z = (4.60) 

At (i = 3 we obtain z = 3/2, which is the expected result for ^He. This 
result was obtained first by Patashinskii and Pokrovskii [SS]- Note that 
the same order in when Tc is approached from above fails to give a 
non-trivial dynamic scaling behavior. Only after taking into account non- 
Gaussian Gaussian fluctuations, corresponding to the next-to-leading order in 
is possible to obtain a non-trivial dynamic exponent. The origin of this 
non-symmetric critical behavior comes from the intrinsic existing asymmetry 
in the dilute Bose gas with respect to the ordered and disordered phases. 
Indeed, in the ordered phase the spectrum has a relativistic-like form, while 
in the disordered phase the non-relativistic behavior dominates the physics. 
The static critical exponents are not affected by this asymmetric behavior 
of the theory, but the critical dynamics properties are. Note that our value 
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of the dynamic exponent is independent of A^, i.e., z = {Q — d)/2. This 
exponent agrees with model F critical dynamics only at d = 3. There the 
dynamic exponent is given exactly hj z = d/2 [70]. At d = 3 the value of z is 
indeed expected to be 2; = 3/2 [57]. However, a word of caution is necessary 
here. Our calculation of the dynamic exponent was made assuming that finite 
temperature dynamics can be derived out of a quantum Hamiltonian using 
equilibrium statistical mechanics. This is not quite right. Critical dynamics 
cannot be derived from a Hamiltonian and equilibrium statistical mechanics 
analysis [70] . 

The full classical theory, including dynamical effects, should account for 
the frequency dependence of fl{iun,p)- The full expression for fl{iun,p) in 
the classical approximation and 2 < d < 4 is [53] 



H(iw„, p) = -Aam Tp 



where 



2miuJn 
p2 



d-3 



+ e-^^^'^-^) ( 1 + 



2miu}n 
p2 ^ 



d-3 



(4.6 



Ad = 22-'^7r-'^/2g*-(rf-2)/2p^^/2 - l)r(3 - d). 



(4.62) 



The derivation of Eq. ( [4.61 ) is made in Appendix [E] As a;„ — )■ Eq. (4.61) 
reduces correctly to 11(0, p) = no(p) = —CidTrn?p'''~^ and we recover Eq. 



(4.57) for the effective interaction. 



For d = 3 Eq. (4.61) becomes 
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fl{iujn, p) 



-i— In I 



2miUn — p 



(4.63) 



By making the replacement iun E{p) + i6 and substituting in Eq. (4.53) 



we obtain after taking the real part of the effective interaction the following 
expression for the excitation spectrum 



E{p) 



PL 
2m 



1 



AmnNg 



P 



,2+ (^TrgV^i^ 



2m_E(p)+p2 



1/2 



(4.64) 



2mE{p)—p^ 

4.2 A Depletion of the condensate 

At T = the depletion of the condensate is more easily obtained from the 
formula 



where 



(4.65) 



= lim 



du 



cfip 



(4.66) 



;To+J_^27tJ {2ny 
Note that for T = the polarization bubble vanishes and we can set t{iu,p) 



g. Explicit evaluation of the w-integral in Eq. (4.66) yields 



d'^p 
(27r)° 



L+^N9 1 
2 



2E(p) 



(4.67) 



where E{p) = ^yp^/Am"^ + nNgp'^ /m. The integral in Eq. (4.67) can be 
easily evaluated using dimensional regularization [H], so that the factor 1/2 
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between brackets does not contribute. Thus, we can rewrite Eq. (4.67) as 



2 V 2 

where = AmnNg and 



') = Ua + ^B], (4.68) 



A= I , ^ :, (4.69) 

' (27r)^^p2 + M2' ^ ' 



and 



(27r)'^p^p2 + jVf2 
The integrals A and -B are related to the integral 



B^l^^^. (4.70) 



by 



■^-''S?^'/.(<'+l). (4.72) 



and 



The integral can be evaluated with the usually tricks already employed in 
Chapter |2] We have evaluated it in Appendix [Bj The result is given in Eq. 



(B.13). Thus, 
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2, {mnNg/Tiy/'^ 



From the relation V{z + 1) = zViyz), we have 



2 \ 2 2 \ 2 \ 2 



(4.74) 



r ( 1 - = r (-i) , (4.75) 



2 2 \ 2 



such that 



2, {mnNg/7iY/^{d-2)^ f d\ ^ f d - 1 



It can be shown that for 2 < d < 4 the couphng g is given in terms of the 
s-wave scattering length by 



(4.78) 



^ T{d/2-l)m' 

For d = 3 this yields the well-known formula g = AiTa/m. Setting d = 3 
finally yields 



8 Nna^\ 

no = n\ (4.79) 

which for iV = 1 is the usual Bogoliubiv's formula for the depletion of the 
condensate 
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4.2.5 The superfluid density 

In order to calculate the superfluid density, we have to perform a Galilei 
boost in the system pH]. A normal fluid is insensitive to a Galilei boost 
but a superfluid is not. This fact imposes strong constraints on the form 
of the excitation spectrum of a superfluid [17] . For instance, it implies that 
an ideal Bose gas in its condensed phase is not a superfluid, although we 
sometimes speak in this case of "ideal superfluidity" , since the calculation of 
the superfluid density using the standard deflnition implies in this case that 
it is identical to the condensate density. This is absolutely not the case in 
general. Indeed, superfluidity may occur even when no Bose condensation 
is possible, as for example in the case of two-dimensional interacting Bose 
systems at flnite temperature, where the Hohenberg-Mermin-Wagner [T9|[20] 
theorem forbids the appearance of a Bose condensate. In this situation there 
is a phase transition to a superfluid state without spontaneous symmetry 
breaking, which is the Kosterlitz-Thouless phase transition [21] discussed in 
Chap. |3} Another example comes from three-dimensional superfluids. It 
is easy to show that for a homogeneous superfluid at zero temperature the 
superfluid density corresponds to the whole density [69] , while the condensate 
is depleted due to interaction effects [ITJ SH] • Experimentally it is estimated 
that for ^He in three dimensions and at zero temperature the condensate 
fraction is about 13% [SH], although the whole system is superfluid. 
A Galilei boost of momentum ko changes the Lagrangian to 

^ko = ^ + |^|&P + ko-j, (4.80) 
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where 



j = (b*Vb-bVb* 

2m 



(4.81) 



is the current density. 

The free-energy density in the presence of the Gahlei boost is given by 
the functional integral 



In^ Vb*Vbe'-^o'^^-^'^''^'^o 
The superfluid density is defined by [58] 



(4.82) 



Ps = rn? lim 



kj^o (9kQ 



It is now easy to show that 



(4.83) 



P. = m(|6r) 



-/ ,lr j ,f rC(r,r) 



(4.84) 



where 



C(r,r) = (j(r,r) ^(0,0)) - (j(r,r)) ■ (j(0,0)) (4.85) 



is the connected current correlation function and the factor l/d'm Eq. (4.84) 
emerges due to rotational invariance. Note that is the total density of 

the fluid and sometimes we define p = m(|6p) as the total fluid density. 
As an example, let us calculate the superfluid density for the Lagrangian 



(4.80) in a regime where the non-Gaussian fluctuations are small using a 
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Bogoliubov transformation [18] for the Lagrangian (4.80). We will not con- 
sider the effects of vortices, which is certainly important when perfoming 
calculations using directly the phase of the order parameter [21 [21] . 
In the Bogoliubov approximation the superfluid density is given by 



rp oo 



E / [G^K,k) - |F(a;„,k)p] , (4.86) 

a=-oo ^ ' 

where we the Green functions are given by the zero temperature counterpart 



of the Green functions (4.47) and (4.48), which we write as 



k) = (4.87) 



and 



^K, k) = ^ .f", (4. 



2 I l^Y _L _2P b-2 
^ \2m) ^ 



In the above equation 6o = {b) and we are approximately setting ~ l&oP, 
which is true up to higher order effects. 

The above integral and Matsubara sum can be calculated exactly if we 
assume a low-energy approximation to the Bogoliubov spectrum, i.e., 

E(k) ^ c|k|, (4.89) 

where (? = gp/rn^. Thus, by rescaling the momenta to remove the c factor, 
the expression for the normal fluid density (note that ps = p — Pn) becomes 
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n=— oo " \ ' " 



2T y 



n= 



c<^+2d ^f-^_ 7 (27r)'^ (a;2 + ^2)2 

which can be rewritten as 



_ T r d'^k T ^ r d^k 

_ V [J^^^ (491) 

n=— cxD >' ^ ' ^ " ' 

It is convenient in this case to perform the integrations first, and only 
afterwards evaluate the Matsubara sums. The first and second integrals in 
Pn can be related to the integrals Ji and I2 evaluated in the Appendix [B| 
where u"^ plays the role of M^. We only need to notice that 



^dh, 1,2 

■-27Td Ia{d + 2), (4.92) 



{2ttY (F + M2)° 

where a = 1,2. Note that we have to replace c? — )■ c? + 2 in the integrals Ii 
and I2 of Appendix [Bj Thus, 



d'^k k2 TT'^/^T'^ f d\ , , , 

-r 1 - - |n|''. (4.93) 



(27r)'^u;2 +k2 d \ 2_ 

Now the sum can be easily performed by using the definition of the zeta 
function (see Appendix [C]) , so that we obtain. 
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oo „ 



d 



V 



(l-0CM)- (4.94) 



The remaining integral along with the Matsubara sum in p„ is easily eval- 
uated using a similar procedure. Therefore, after some simplifications, the 
final result for the superfluid density is 



where A is an ultraviolet cutoff which is of the order of the inverse scattering 
length, A ~ . The term linear in the temperature arises also in a mean- 
field calculation of the stiffness in the classical XY model [60] ■ The term 
proportional to T'^'^^ is a quantum correction due to phonon excitations of 
the Bose liquid [17] . The Landau formula usually only gives the contribution 
proportional to T'^'^^. 
For d = 2 we have 



2m'^+2 r j^drp 



Ps 



^ {gpY+'^/^d [{ATTy/^dT{d/2) 
(d + l)7r'^/'r (l-^) (:{-d)T'' 



(4.95) 



+ 3C(3)r 



i3 



(4.96) 



P 



2n{gpy [ 4 



while ior d = 3 we obtain 




(4.97) 
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It should be noted that for d = 2 the couphng g itself is density dependent, 
being given by |61] 



47r 

(4-98) 



4m 



where 7 is the Euler constant and the naive diluteness condition pa'^/ m ^ 1 
has to be modified for d = 2 to lnln[m/(pa^)] ^ 1 [62j. This is in contrast 
with the d = 3 case where g = ATra/m, with no dependence on the the total 
density. 
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Chapter 5 
Mott insulators 



5.1 The Hubbard model 

The simplest electronic model of a Mott insulator is provided by the so called 
Hubbard model: 

^ ^-t^Yl •fi'^-f^" - I^^^riia + UY^ rii^riii, (5.1) 

a i,a i 

where fia- is a fermion in the lattice, rii^ = fj^fia, and a =t)4- the above 
Hamiltonian t,U > 0, and is the chemical potential. In the kinetic term 
the lattice sum runs only nearest neighbors only. The Hubbard Hamiltonian 
is exactly solvable in the limits U — and t — 0. The non-interacting hmit 
corresponds to band theory in the tight-binding approximation, while the 
t — limit corresponds to the atomic limit. Although in the atomic limit 
the theory is interacting, it is easily diagonalizable, since 
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Ht=o 



■11 



(5.2) 



where 



(5.3) 



Thus, 



[hi,nia] = 0, 



(5.4) 



and therefore the eigenstates of also eigenstates of hi. In this case 

we can simply omit the lattice sites of any calculation, since the sites are 
decoupled. The eigenstates of the h are |0), If), \ i), and | '\D, correspond- 
ing to empty, singly occupied (with either up or down spins), and doubly 
occupied sites, respectively. The corresponding eigenenergies are Eq = 0, 
= El = — /i, and 62 = U ~ 2ji. 

Let us consider the example of a half-filled band, i.e., the total number 
of fermions in the system equals the number of lattice sites L. In this case 
it can be shown that for a bipartite lattic^ /i = U/2 exactly. This result 
is straightforwardly checked in the atomic limit. To see that this result is 
also valid when t 7^ 0, we just perform a particle-hole (ph) transformation 



fi, ^ e^Q-^7iL, ft ^ e^Q-^7i<x, where Q = (7r,...,7r). When /i = U/2 the 



Hamiltonian is invariant under this particle-hole transformation. Note that 
the factor e**^ "'^' is necessary in order to maintain the sign of the hopping 

-'^This means that the lattice can be viewed as being made of two interpenetrating 
sublattices. One simple example is the cubic lattice 
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term, because fermions anticommute. Indeed, we have that g^Q CR-i-r'-j) = 
since Rj and Rj are nearest neighbor sites. Thus, the ph-transformation 
yields in general the Hamiltonian 



H' = U-2^i-tJ2Yl flf^^ + {^^-U)J2n^. + uJ2 ^^t^4 (5.5) 

(ij> a a i 



Note that the Hamiltonians H from Eq. (5.1) and the one given above 



coincide for = f//2. Now let / and /' be the free energy densities of the 



Hamiltonians H and H' given in Eqs. (5.1) and (5.5), respectively. We have 



the particle density is given as usual by the thermodynamical relation 



n 



djj,'' 



(5.6) 



while from /' we obtain. 



n = — 



(5.7) 



Since ii = U/2 implies H = H' , we also have that f = f for this value of 



the chemical potential. Therefore, for = U/2 the RHS of both Eqs. (5.6) 



and (5.7) are the same, so that we obtain n = 2 — n, which yields n = 1. 



At half-filling the Hubbard Hamiltonian can be rewritten as 



{id) * 



(5.8) 



where 
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Si = ^ XI fl^^pf^^^ (5-9) 



2 



with cr = (cTi, (72, cTs), cTj being the Pauh matrices. For U ^ t doubly occupied 
sites are strongly suppressed and second-order perturbation theory yields the 
effective Hamiltonian [3] 



4^2 

H = — S, ■ S„ (5.10) 



subjected to the local constraint 



J^rii. = 1. (5.11) 



Note the subtlety here. The half-filling condition demands that the particle 
density 



n 



yX(n,.) = l, (5.12) 



which is easily enforced when /i = U/2. This is a global constraint, that 
simply demands the average site occupation be the unity. However, when 
U ^ t the average constraint becomes a local one given by the operator 



equation (5.11), i.e., double occupation is strictly forbidden. 



The effective Hamiltonian (5.10) is the one of a Heisenberg antiferromag- 
net. It is rotational invariant in spin space, i.e., it has an SU{2) symmetry. 
This means that the total spin operator 
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5^ Si, (5.13) 



commutes with the Hamiltonian. However, the ground state of the Heisen- 
berg Hamiltonian breaks this symmetry. The most favorable state at zero 
temperature corresponds to alternating spins in the lattice. This state is pic- 
torially shown for a square lattice in Fig. 2.1. The antiferromagnetic state 
shown in the figure is actually a mean-field state for the Heisenberg model, 
the so called Neel state. This state is also a mean-field state of the Hubbard 



model at half-filling. Incidentally, the total spin operator (5.13) obviously 
commutes with the Hubbard Hamiltonian, showing that the Hubbard model 
is SU (2) symmetric, as expected physically. 

Let us perform the mean-field theory for the Hubbard model in a d- 
dimensional cubic lattice explicitly. In order to do this, we introduce an 
auxiliary field via a Hubbard-Stratonovich transformation: 

^ = -tEE/U.-f^E--S. + ^$^m,l (5.14) 
We are looking for a mean-field state with a staggered magnetic moment 

mi = e^'^-^'m, (5.15) 

where the vector m is uniform. Rotational invariance allows us to fix a 
direction for m. We will choose the quantization axis to be along the z- 
direction. Thus, m = me^;. Let us define 
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f — ^ — f — ^ 

/<■ 

)^ Jif 

Figure 5.1: Mean-field ground state for a Heisenberg antiferromagnet, tlie so 
called Neel state. 
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fia — 



Cia, i e B 
The mean-field Hamiltonian can be written as 



^ 3C/L 



k,(T 



where 



Ck(T 



Ck(T 



and 



ko- 



aUm 



2 



with the tight-binding dispersion 



£k = — 2t cos /C(j. 



a=l 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



The mean-field Hamiltonian is easily diagonalized and leads to the energy 
spectrum 



(5.21) 



and we see that the electronic spectrum is gapped. Thus, the mean-field 
ground state energy per site is 



^o = -|E^k^ + f-^^ (5.22) 

where the prime on the sum over k is to denote that we are summing over 
the upper half of the Brillouin zone. Let us speciahze to two dimensions. 
By minimizing the above equation with respect to m, we obtain the gap 
equation, 

3 f dkrj. f dky 1 . . 

In two dimensions it is not a too bad approximation, at least for more qual- 
itative purposes, to use a square density of states to evaluate the above 
momentum integral, 

where W = 4t is the bandwidth and 6{x) is the Heaviside function. This 
yields the magnetization, 

2W e-w 

^=~iT~. (5.25) 

u I — e u 

Thus, mean-field theory predicts that at half-filling the Hubbard model is 
an antiferromagnetic insulator for all f/ > 0. This is not quite accurate, of 
course. For a small enough U the system should be a metal. Indeed, the 
Hubbard model undergoes a metal insulator transition for a critical value of 
the interaction typically of the order of the bandwidth W ^64j- The metal- 
insulator transition in the Hubbard model started to be better understood 
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with the use of the so called dynamical mean-field [65] , which is based on the 
large dimension limit of the Hubbard model [66j. For a thorough review on 
the subject, see Ref. [65] . 



5.2 The Bose-Hubbard model 

The Hamiltonian of the so called Bose-Hubbard model [HZ] is given by 

where U, J > 0, hi = h\hi and /i is the chemical potential. The summations 
are over the sites of a cubic lattice and the symbol means a sum over 
nearest neighbors. The operators hi obey the usual bosonic commutation 
relations, i.e., = 6ij und = [^J,^j] = 0. 

The aim of this tutorial is to provide an introduction to the theory of the 
Bose-Hubbard model, which in the last years gained considerable experimen- 
tal relevance in the context of Bose-Einstein condensation (BEC) [68j . 

What kind of phases we expect for the above model? Firstly, let us note 



that the Hamiltonian (5.26) is simply a lattice version of the interacting Bose 
gas Hamiltonian. Indeed, the hopping term is just a lattice derivative. Thus, 
for small enough U we expect to obtain a superfiuid featuring the well known 
Bogoliubov spectrum. To see this, just write 

k = ho + 6K h\ = hl + 5hl (5.27) 
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where bo represents the condensate and minimize the Hamiltonian, and 5bi 
are small fluctuations around the condensate. By keeping just the quadratic 
fluctuations, it is easy to see that the Hamiltonian can be approximately 
written as H — LEq + 5H {L is the number of lattice sites), with 



Eo = - {2dJ + n + -j no+-nl, 



where we have performed a Fourier transformation in the lattice and 



(5.28) 



(5.29) 



Sbl 5b. 



5b^ 

6P „ 



(5.30) 



Ml 



Ubl 



(5.31) 



Ek- fx-U/2 + 2Uno 

U{b*f £k - /X - t^/2 + 2Uno 

with £k = ~2J^^^^ cos and Uq = |&oP- Due to the minimization condi- 



tion we have 



fj, — —2d J — — + UriQ. 



(5.32) 



The above result follows easily by demanding that the linear terms in both 
56j and vanish. Alternatively it may be derived by simply minimizing 
with respect to tt-q. The energy spectrum can be obtained by simply solving 
the Heisenberg equations of motion. To this end we need the equations 



i9t56k = [(^&k, ?>il\ idtSbl^ = [56Lk, ^H]. (5.33) 
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After straightforward evaluation of the commutators, we can rewrite the two 
equations above as a single matrix equation: 



tasdt^i, = Mk^k, (5.34) 
where is the third Pauli matrix. The Ansatz 

Mt) = e-'^''*^k(0) (5.35) 



solves Eq. (5.34) provided det(£'k0'3 — Mk) = 0, or 



Ek = ±V(£k - - U/2 + 2f/no)2 - f/%2. (5.36) 



Note that due to the form of the chemical potential (5.32), the above spec- 
trum is gapless, i.e., -Ek=o = 0, as required by superfluidity. 

The above results are valid for U small, i.e., U J. So, what happens 
now in the opposite limit, when U ^ Jl This is a subtle question. In the 
continuum limit exact arguments [011 EO] involving the Ward identities show 
that at zero temperature the superfluid density is identical to the particle 
density, and this for for all values of U . This means that at zero temperature 
the whole system is in a superfluid state, although not every particle is con- 
densed, since the condensate is depleted due to interaction effects. Does this 
exact result also holds in the lattice? The answer is: it depends on whether 
the particle density (n) is integer or not! For noninteger particle density, ap- 
proaching the strong coupling limit from the weak coupling one by varying 
J /U essentially does not change the superfluid characteristics of the system. 

92 



Thus, in this situation the system is still a superfiuid for U ^ J. However, 
for (n) = n eN the situation is different. In this case the bosons will localize 
for large enough U and the system will become an insulator, whose ground 
state has n particles per site. In order to better understand how it actually 
works, let us consider the eigenstates of the number operator (we omit the 
site index for simplicity) 



\n) = 4t(6^)"|0), (5.37) 



and the coherent state 



1^) = e-l^l /2exp(^6t)|o) 

oo „ 

= e-l^l^/^5:^|n). (5.38) 



n=0 



Recall that for a coherent state b\z) = z\z). Let us consider for simplicity 
a two-site problem. For this case let us assume that the system is in the 
coherent state 

\^) = \zi,Z2). (5.39) 
The expectation value of the Hamiltonian in this state is the energy 

E{zi, Z2) = imi^) = -J{ZIZ2 + Zlzi) + 



i=l,2 



I 1 2 ^ I 1 2 / 1 1 2 -1 \ 

-A*Pi| + —\Zi\ [\Zi\ - 1) 



2 



(5.40) 

In terms of Zi = ^/^e**^', where fii is the mean particle number at the site i. 
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Eq. (5.40) becomes 



1=1,2 



(5.41) 



The above energy is similar to the classical Hamiltonian of a Josephson junc- 



tion [72]. In order to explore this similarity further, we rewrite Eq. (5.41) 
as 



E(An,A^) = -Jv/iV2 - (An)2cosA^ + ^(Afi)2 + E^, (5.42) 



where = + ?7,2 is the total number of particles of the system, An 
ni - n2, Ayj = '^\ - <y52, and 



AT 



^+2r+4^- 



(5.43) 



The energy (5.42) corresponds precisely to the Hamiltonian describing a 



two-level Bose-Einstein condensate via a Josephson junction as discussed 
by Leggett [73] . In this case the variable An plays the role of the momentum 
conjugated to Aip. Therefore, the Josephson current is given by 



(5.44) 



The Josephson effect implies superfluidity. Semiclassically, since An and Aip 
are canonically conjugated, we have the uncertainty relation [72] 



AnAif ~ 1. 



(5.45) 
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The so called Josephson "phase-voltage" relation is here generalized to 



U J An 

Recall that in the case of superconductors we have dfAyD = 2eV/h, where 
the voltage V the same as the difference of chemical potential across the 
junction. In the above equation the role of 2eV/h is played by UAn/2. The 
second term is absent in the Josephson relation. This term appears in the 
context of Josephson junctions in Bose- Einstein condensates |71j . 



If U ^ J, the second term in Eq. (5.42) will constraint ni ^ n2 and a 
commensurate situation ni = ^2 = n = 1, 2, 3, . . . will be favored. It is then 
clear that the current will be zero and we have an insulator. 

Now we will solve the full model approximately using a Green function 
method. The aim is to compute the Green function 

G,,{t) = -z{T[b,{t)b]{Q)]), (5.47) 

where T[k{h)b]{t2)] = ^(ti - t2)&i(^i)^i(^2) + ^(^2 - ^l)&](^2)^^(^l) and e{t) is 
the Heaviside function. The Green function can be calculated exactly in both 
limit cases J = and U = 0. We will consider a solution that corresponds to 
a perturbation expansion in J/U, i.e., around the limit case J = 0. Thus, we 
have to compute the exact Green function in this limit in order to proceed. 
The J = limit is actually a single-site system, since the Hamiltonian 
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can be written as a sum of single-site Hamiltonians: 



where 



Hj=o = J2k (5.48) 



hi = -fiUi + ^hi{hi - 1). (5.49) 



In this case it is enough to compute the Green function for the single-site 
Hamiltonian hi and the site index can even be omitted. Since [h, h] = 0, the 
eigenstates \n) oi h are the exact eigenstates of h with eigenvalues 



E^^ -Hn+^n{n- 1). (5.50) 



We want to compute 

= -immb\o))+e{-t){b\o)bm. (5.51) 

Let us assume a ground state with n e N particles per site. Then we have 

{b{tp{0)) = {n\b{tp{0)\n) 

= Vn + l{n\e'H(0)e-''^*\n + 1) 

= (n + l)e*(^"-^"+i)* (5.52) 

Similarly we find 

{b{Oyb{t)) = ne^(-^"-i-^")*. (5.53) 
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Therefore, 



g{t) = -i[e{t){n + i)e^(^n-^n+i)* + ^(-t)ne'(^"-^-^")*], (5.54) 



such that the Fourier transformation 



oo 

itot 



g{u) = / dte^'git) (5.55) 



reads 



n + 1 n 



U + En — En+l + iS U + En-l — En — i6 

n + 1 n 



(5.56) 



u + fi — Un + i6 u + fi — U{n — 1) — i6' 
where 6-^0^. Note that in order to perform the Fourier transformation 



of (5.54) a convergence factor e""^* was used for the first term, while for the 
second term a convergence factor e*^* was needed. 

The approximation we are going to employ is equivalent to a well known 
mean-field theory approach to solve the Bose-Hubbard model [6]. Instead 
performing a mean-field theory in the Hamiltonian, we will compute the 
Green function approximately using an expansion on the hopping. Thus, the 



unperturbed Green function will be given by Eq. (5.56). In order to better 
motivate the method of solution, let us show how it can be used to solve 
the exactly solvable limit U = 0. For this particular case the Hamiltonian is 
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easily diagonalized by a Fourier transformation: 



^(£k - /^)nk, 



(5.57) 



k 



where 



fej^fek- The exact Green function is obviously given by 



1 



(5.58) 



uj + II — e^ + iS 



Let us derive the above Green function in a less direct way, namely, via 
the hopping expansion. It will be a more or less complicate way to derive a 
straightforward result, but it will serve the purpose of illustrating the strategy 
to solve the Bose-Hubbard model approximately. 

The U = Hamiltonian can be decomposed in the following way: 



Hu=o — Ho + Hi, 



(5.59) 



where 




(5.60) 



is the unperturbed Hamiltonian, and 




(5.61) 



is the perturbation. The unperturbed Green function is given by 



1 



(5.62) 



CO + II + iS 
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The exact Green function can be obtained by performing the infinite sum of 
diagrams shown in Fig. 2.2. The continuum hne represents the local Green 



function (5.62) at a lattice site i, while the dashed line represents a hopping 
process between two neighboring sites. The exact Green function gives the 
response of the system to a boson propagating from a site i at some time t 
to another site j at an earlier time t'. The diagrams of the figure illustrate 
all the possible processes for a quadratic Hamiltonian. In this case only tree 
diagrams contribute are nonzero, since U = 0. The perturbation expansion 
reads simply 

Gij{u) = go{^^)Sij + Qo{u;)JijQo{u) + Qo{^^)^JiiQo{^^)JijQo{uj) + ..., (5.63) 

where Jij = — J if are nearest neighbors and zero otherwise. This has 
the structure of a geometric series and can be rewritten as 

Gij{uj) = go{uj)5ij + go{uj)'^Jii[5ijgo{u) + go{uj)Jijgo{uj) + . . .] 

I 

= go{uj)6,,+go{io)J2j^iGiAuj)- (5.64) 
The Fourier representations 

J., = ^E^^''^'''"''^^^'^' (5.65) 

k 



G^A^) = ^ J]e^'^-(^-^^)G(k,a;), (5.66) 

k 

5,, = ^5^e^^-(^-^^), (5.67) 
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Figure 5.2: Hopping expansion for the Green function. 



allow us to rewrite Eq. (5.64) in the simpler form: 



(5.68) 



which can easily be solved to obtain once more the exact Green function for 
U = 0: 



G{k,uj) 



1 



1 

u + — Ek + i6 



(5.69) 



Unfortunately, when U ^ the series of diagrams given in Fig. 5.2 do 
not lead to the exact Green function, since loop diagrams containing higher 
order local Green functions (higher order cumulants) are missing. The latter 



vanish when U = 0. Nevertheless, the diagrams of Fig. 5.2 still give a good 



approximation, especially in three dimensions. Thus, we can approximate 



the Green function for the Bose-Hubbard model using Eq. (5.64) with go{uj 



replaced by Giuj) given in Eq. (5.56), i.e.. 



G{k,u) 



(5.70) 
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The energy spectrum is given by the poles of the above Green function: 

E±(k) = -i_i+^[s^ + (2n - 1)U] ± ^^4 + 2(2n + l)t/£k + t/2. (5.71) 

For large enough U there is an energy gap A between the + and — branches 
of the spectrum, which is given by 

= V(2rfJ)2 - 4rf(2n + 1)UJ + C/2. (5.72) 

The presence of the energy gap indicates that the system is an insulator. As 
U gets smaller, it will eventually attains a critical value Uc below which the 

system becomes a superfluid. This critical value of U is found by demanding 
that the gap vanishes for U — Uc- The condition A = yields 



= 2dJ 



2n + l±2Vn(n + l)] . (5.73) 



Note that both solutions are positive. In order to know what is the right 
one we have to plot the phase diagram and study it more carefully. The 
transition from the insulating phase to the superfluid phase occurs when the 
bosons condense. This happens when the Green function is singular for w = 
and k = 0, since the bosons condense at k = 0. Thus, the phase diagram is 
given by the equation 

g{Q) = — , (5.74) 

£k=0 
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or 



n 



n + 1 



U 



^i/U + l-n fi/U-n 2dJ' 



(5.75) 



In Fig. |5.3| we plot the phase diagram for d = 3 and n = 1,2,3,4. It 
features the so called Mott lobes [HZl- The phase inside the Mott lobes is 
an insulating one. Outside them we have a superfluid phase. Each Mott 
lobe corresponds to a Mott-insulating phase with n particles per site. Thus, 
the largest lobe, corresponding to < fi/U < 1 has n = 1. The next one, 
in the interval 1 < yu/f/ < 2 has n = 2, and so on. The tips of the lobes 
correspond to the points where the upper and lower bands of the spectrum 
meet for k = 0, thus closing the gap. The coordinates of the tips can be 



easily obtained by extremizing Eq. (5.75) with respect to fi/U. The tip of a 
Mott lobe corresponds to the maximum value of J/U for a given value of n. 
Extremization gives the results: 



v/n(nTl) - 1, 



(5.76) 



f/c = 2dJ 



2n + 1 + 2^/n{n + 1) 



(5.77) 



and we see that Uc corresponds to the solution of the equation A = 0. 

Let us compare the above results with recent highly precise Monte Carlo 
(MC) simulations for the n = 1 case [71]. There it is obtained that J/Uc ~ 



0.03408. On the other hand, our formula (5.77) for d = 3 and n = 1 gives 



J 



Uc 6(3 + 2^2) 



0.0286. 



(5.78) 
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Figure 5.4: Phase diagram for n = 1. The Black circles are Monte Carlo 
points from Ref. The continuous line is the mean-field result. 

Near the tip of the Mott lobe, which corresponds to the quantum critical 
regime, our approximation is very bad. This is to be expected, since our 
approach is equivalent to mean-field theory. In Fig. |5.4| we compare our 
result for n = 1 with the MC phase diagram of Ref. [71]. Indeed, the MC 
points agree with the mean-field curve only far away of the critical point, i.e., 
for small J/U. 

A better approximation, where higher order cumulants are included, so 
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that also loop diagrams appear in the hopping expansion, allows to go beyond 
mean-field theory and approach better the MC result. This was done recently 
in a remarkable paper by dos Santos and Pelster [TSj, where a quantum 
Landau-Ginzburg formalism including higher order cumulants was developed. 
Their result closely agrees with the MC simulations. Furthermore, in contrast 
with the mean- field approach discussed here, their analysis shows how the 
behavior of the phase diagram changes with the dimensionality. Their phase 



diagram is shown for = 3 in Fig. 5.5 



5.3 The Heisenberg antiferromagnet 
5.3.1 Antiferromagnetic spin waves 

We have seen in Chap. [2] that the spin waves in a ferromagnet have a dis- 
persion u ^ k"^. For an antiferromagnetic system the situation is different. 



due to sublattice magnetization (see Fig. 5.1). Thus, in order to derive the 
spin-wave excitation of a Heisenberg antiferromagnet we have to consider the 
lattice sites explicitly. We will perform the calculations in the semi-classical 
limit. The Heisenberg model for an antiferromagnet will be written as 

= JS^^ni ■ iij, (5.79) 

(id) 

where = 1. We can easily write down the corresponding lattice Landau- 
Lifshitz equation if we rewrite the above Hamiltonian as 
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Figure 5.5: Phase diagram of the Bose-Hubbard model in three dimensions 
beyond the mean-field approximation discussed in the text. The result is 
from Ref. [75]. The solid green line corresponds to the effective action 
approach employed by dos Santos and Pelster |75j. The dashed purple line 
corresponds to a variational approach used by the same authors. Both are 
compared with the MC results from Ref. [71] (red triangles), mean- field 
theory (dotted dashed blue line), and third-order strong-coupling expansion 
(dotted black line) [75] . 
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H = JS^^ni- Bi, (5.80) 

i 

where the effective local magnetic field is given by 

Bi = Y,n„ (5.81) 

j 

where the prime in the sum reminds us that the sum is over the sites which 
are nearest neighbors of i (that is why the effective field is i-dependent). For 
example, in one dimension we have 

Bi = ni_i + rii+i. (5.82) 
Therefore, the LL equation is given for the lattice model as 

dtUi = JS^Ui X Bi). (5.83) 

For simplicity, we will derive the spin- wave spectrum in one dimension. 
The generalization to higher dimensions is straightforward. 

In one dimension the staggered magnetization occurs between even and 
odd lattice sites, such that we have 

n2m = ^^063 + Sn^^, (5.84) 

and 

n2m+i = -^1063 + (5-85) 
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where 6nj- is a small transverse fluctuation that will be used to linearize the 
LL equation. Thus, the calculation is similar to the one for ferromagnets in 
Chap. [2| except that here we have to consider a staggered magnetization 
and work on the lattice. Thus, we have the coupled equations, 



(5.86) 



dt^2m+l — JS'^[Tn-2m+l X (n2m + n2(m+l) )] ■ 



(5.87) 



Substituting Eqs. (5.84) and (5.85) in the above equations and neglecting 



terms of order higher than two in the fluctuations, we obtain. 



(5.88) 



and 



''2m.+l)i 



(5.89) 



where we have defined dri^ = 6n^ + i6n'j. 



We will solve the Eqs. (5.88) and (5.89) via the Ansatze: 



^ntm+i = ve 



i[{2m+l)k-Ljt] 



(5.90) 



Thus, Eqs. (5.88) and (5.89) can be put in matrix form 
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u 












V 







U 




V 



(5.91) 



—1 —cos A; 
cosk 1 

where uq = 2noS'^J. The frequency cu is therefore determined by the eigen- 
value equation, 



det 



— {uj + ujq) —ujo cos k 

OUq cos k UJq — UJ 

which yields the spin- wave spectrum: 



0, 



(5.92) 



(k) — oJoV 1 — cos^ k — ojo \ sin A;|. 



(5.93) 



The generalization to higher dimensions is straightforward, leading to the 
result. 



a;(k) = ojQ 



(P — j cos ka 



,a=l 



In the long wavelength limit, the above equation becomes 

uj{k) \/d uJo\k\. 



(5.94) 



(5.95) 



Therefore, in contrast with a ferromagnet, which has a quadratic dispersion, 
the antiferromagnet has a linear dispersion. 
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5.3.2 The quantum 0{n) non-linear a model 

In Chap. [2] we have studied the classical 0{n) non-linear a model in de- 
tail. In that Chapter the model was classical because it was being used to 
describe the thermal fluctuations in a ferromagnet. However, in the context 
of an antiferromagnet, the classical 0{n) non-linear a model is used to de- 
scribe the quantum dynamics of an antiferromagnet at T = 0. In this case, 
the temperature of the classical model studied in Chap. [2] is replaced by a 
coupling constant g and the dimension d should be interpreted as D + 1, 
such that d corresponds to the number of dimensions of spacetime with time 
being imaginary. In this context D is the dimension of space. 

The action of the quantum 0{n) non-linear a model is given by 



5 = 1 dr y rf^x[(9,n)2 + (diu) ■ {dm) + zA(n2 - 1)]. 



(5.96) 



The zero temperature case follows directly from the results for the clas- 
sical model, with T replaced by g, and d being interpreted as the dimension 
of spacetime. So, let us consider the so called quantum critical finite tem- 
perature case. In this regime, we have that g = gc, but the temperature is 
not zero, so that the mass gap m is no n- vanishing. In fact, in this case the 
temperature is the only energy scale available (up to the energy cutoff) and 
we actually expect that = aT^, where a is a number to be determined. 
For g = gc and T > we still have s = 0, so that the relevant gap equation 



is (2.113) (with T replaced by g), which at finite temperature becomes 



110 



^ I d^k 1 _ 1 _ /" d^^^P 1 



where we have used Eq. (2.111) with T replaced by (7. We stress that the 
temperature T appearing in the above equation has nothing to do with the 
temperature T in Chap. \^ as in the present case T arises from the integration 
over r G (0, /3 = 1/T) in the action. 



We can perform the Matsubara sum straightforwardly (see Appendix D ) 
to obtain 



rf^+V I f d^k 1 If d^k 



(5.98) 

Note that by considering the T = limit, we have from the above equation 
that 



fiD+i^ 1 r d^k ^ 

dpi I d k I 



(27r)^+ip2 J (27r)«2|kr 
The integrals can be easily evaluated for D = 2 and we obtain 



m = -2rin(l-e-'"/'^), (5.100) 



whose solution is 



m(r) = rin I 1 . (5.101) 
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The above result has been obtained before by Chubukov et al. [77]. 



5.3.3 The CP^ ^ model 

Many two-dimensional Mott insulators feature competing orders. One exam- 
ple is the competition between a Neel state and a valence-bond solid (VBS) 



state [79J, as illustrated in Fig. 5.6 In a Landau-Ginzburg- Wilson (LGW) 
framework, competing orders usually feature a first-order phase transition. 
This is because the order parameters involved are the most fundamental fields 
in a LGW type of theory. Furthermore, within this so called LGW paradigm 
[78l [8T| l82] second-order phase transitions are characterized by a very small 
anomalous dimension of the order parameter. In quantum phase transitions, 
on the other hand, there are examples of order parameters which are them- 
selves made of more elementary building blocks. This is precisely the case of 



the Neel- VBS transition illustrated in Fig. 5.6 In this case both order fields, 
the staggered magnetization orientation field n and the valence-bond order 
field ipYBSi are made of spinous, which in this case are more conveniently 
represented as a bosonic excitation. 

Mott insulators featuring competing ordered states have been matter of 
intensive discussions in recent years [78j. In particular, there has been some 
controversy on the nature of the Neel- VBS phase transition [83] EU [851 ESI 
l87t 188] . The LGW point of view ignores the topological order of the Mott 
insulator. This topological order is essential in the VBS phase, since degen- 
erate VBS ground states are connected by tunneling events characterized by 
a topological charge [72] . This important aspect of Mott insulators is more 
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Figure 5.6: Schematic phase diagram showing a (second-order) quantum 
phase transition between a Neel state and a VBS as a function of a dimen- 
sionless couphng g. The different order parameters are shown. In the Higgs 
phase the spinous condense due to a spontaneous U{1) symmetry breaking. 
This Higgs mechanism produces an antiferromagnetic state. In the confine- 
ment phase the excitations are gapped and the spinous are confined. While 
in the Neel phase the emergent photon is gapped, in the confined phase it is 
the dual of the photon which is gapped. 

easily seen by performing a CP^ map of the unit vector n giving the direction 
of the magnetization into a larger space defined by some set of complex fields, 
i.e.. 



n = zl(TabZb, (5.102) 

where + \z2\'^ = 1 and cr = (cti, (T2, (T3) is a Pauli matrix vector. Since 
= 1, we have a map of the two-dimensional sphere 5*2 into a three- 
dimensional sphere S3 (since + |z2p = + + + /3| = 1, with 
Za = o^a + iPa)- It reprcscuts n as a composite field having the elementary 
constituents Za, the so called spinous. One important aspect of the CP^ map 
is that it introduces a gauge structure in the system. Indeed, the field n 
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written in terms of spinons fields is a gauge invariant object, since the lo- 
cal phase transformation Za{x) — ?■ e''^^^^Za{x) leaves n invariant. This gauge 
symmetry can be better understood by considering the topological charge 

Q = ^ (p dS^e^^xT^ ■ {d^n x ^aIi), (5.103) 

where Q G N. We have encountered the above topological charge before in a 
different context in Chap. |2| Eq. (2.31). There it was a spatial topological 



object that we called a "hedgehog", a magnetic monopole-like excitation. 



The main difference between the topological charge in (5.103) and the one in 



Eq. (2.31) is that the former lives in a (2 + 1) -dimensional spacetime, while 
the latter arises in three-dimensional space. The spacetime hedgehog is often 
called instanton. 

In order to obtain the gauge field from the CP^ representation we need 
the formula 



(^abcCr1,^<ya(sO'!yS = 2i{6 f,s5aiy5 13j - 6 ^pdu-ySas) , (5.104) 

which can be easily derived by using the commutation relation [a", a''] = 
2ieabcO''^ and the completeness relation 

Cr^pO'^S = '^^a5Sl3^ - SapS-yS. (5.105) 

Therefore, we obtain from the CP^ map that the topological charge is given 
by the magnetic flux 
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27r 

where = e^^^xdiyAx and 



Q=lL'P dS^B^, (5.106) 

S2 



A,= '-{z:d,Za-Zad,Z:). (5.107) 

By coupling the above gauge field minimally to the spinous we obtain the 
Lagrangian for the CP^ model 



C=-\{d^-iA^)za\\ (5.108) 
9 



In view of Eq. (5.107) and the form of the above Lagrangian together with 



n = zl^cTabZb along with the corresponding constraint, it is not difficult to 



see that the Lagrangian (5.108) is equivalent to the Lagrangian of the 0(3) 
nonlinear a model. 

The J — Q model: A lattice model for the AF-VBS transition 

Competing orders need competing interactions. For instance, the AF-VBS 
transition can be obtained from the following lattice model due to Sandvik 
ESI: 



(ij) (ijkl) ^ ^ ^ ^ 



(5.109) 



where the sum over ijkl is around a plaquette, as shown in Fig. 5.7 The 



phases of the model are as shown in Fig. 5.6 When Q/ J -^1 the Heisenberg 
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5.. 



#5/ 



Figure 5.7: Schematic representation of the plaquette term in the J — Q 



model, Eq. (5.109). The sum goes over the pairs of spins indicated by the 



red bonds at left and right panels. 

term dominates over the Q-term and a Neel state is favored. For Q/J ^ 
1, on the other hand, the Q-term determines the ground state, favoring a 
crystalline pattern of singlet valence bonds. Recall that 



Pi 



- S,: ■ S 



J' 



(5.110) 



is a projection operator for singlet states. Indeed, we can rewrite the above 
operator as 



J 



■'3/ 



(5.111) 



Note that when the total spin of the (i, j)-bond is zero, Pij yields the unit as 
eigenvalue. For a triplet bond it vanishes. 

There is some discussion in the literature on whether the J — Q model 
really features a quantum critical point or, in other words, on whether it 
undergoes a second-order phase transition [HH EHl [M] . A related issue refers 
to the gauge theory description of the J — Q model in terms of a CP^ model 
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with instanton suppression via a Maxwell term [78l [82], [83l [HU [85]. While 
Kuklov et al. [85j numerically find a weak first-order phase transition for 
the CP^ model in 2 + 1 dimensions with a Maxwell term, a recent paper 
by Sandvik [89] shows strong evidence that a second-order phase transition 
occurs in the J — Q model. However, logarithmic violations of scaling are 
also found, which may cast some doubts on whether CP^-like models can 
really fully provide a description of quantum criticality in the J — Q model. 

The CP^^ model in the large N limit 

The CP^ model can be generalized to include N spinous instead of two. This 
yields the so called CP^~^ model, which can be studied in the large limit. 
In order to facilitate the analysis, we write the CP^^^ model as 



such that the constraint is accounted for by the introduction of a Lagrange 
multiplier field a. Note that g here is a dimensionless coupling. We have also 
introduced a Maxwell term in the Lagrangian. With this Maxwell term the 
model at = 2 is no longer equivalent to the 0(3) nonlinear a model. This 
occurs only in the Cq — > c>o limit. By integrating out N — 2 spinon fields we 
obtain the effective action 




[\id,-tA,)z, 



+ 



(5.112) 
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= {N-2)Ti\n\-(d^-iA^f + ia] 



+ iai\za\' - 1)]} . 



(5.113) 



The saddle-point equations are obtained in a standard way. We will consider 
a saddle-point with = 0, ia = o-q, and Za = (a, where ctq and (a are 
constants. This means that no classical instanton solution is being considered 
here. 

The calculation below is very similar to the one made for the classical 
0{n) nonlinear a model in Chap. [2| so the reader may be interested in going 
back to that Chapter again to recall the main features of the calculation. 

The VBS phase occurs for g > Qc, where Qc is a critical coupling. In this 
case we have ctq 7^ and (a = 0, leading to a gap equation (for large A^) 



For g < gc, on the other hand, we are in the Neel phase and ctq = and 
(a 7^ 0. Thus, we have 



For g = g^we have ctq = and (a = 0, which determines the critical coupling: 
gc= (d- 2){2nY/{NSd), where Sa = 2t:'^/^ /V{d/2) is the surface of a unit 
sphere in d dimensions. 




(5.114) 




(5.115) 
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It is easy to see that in general we should have 



E + - - 1 = NgK--'a, f ,f^, ,\ (5.116) 



Therefore, we obtain for g > gc, 



9c 



-| 2/{d-2) 

(5.117) 



r(d/2)r(2 - d/2) V g 

It is clear that (Jq should be identified with where ^ is the correlation 
length. This leads to the well known large N critical exponent for this quan- 
tity, u — l/{d — 2). The limit d — > 2, on the other hand, yields 

ao = A2exp(^-^^, (5.118) 

and we see that in this case there is a gap for all > 0. 

For g < gc we obtain X]a=i2 iCal^ — 1 ~ 9/9c, which implies that the 
critical exponent of the Neel field order parameter 

(n) = C>abCb (5.119) 

is — 1- Thus, from the hyperscaling relation /3n — u{d — 2 + r}N)/2 it 
follows that the anomalous dimension is 

r]N^d-2 (5.120) 

for large N. 

Now we consider the gauge field fiuctuations to lowest order for g > ge- 



lid 



This correction is obtained from the expansion of 

TT\n[-{d^-tA^f + ao] (5.121) 

up to quadratic order in the gauge field. In order to do this, we first note 
that the covariant differential operator appearing in the inside the logarithm 
can expanded as 

- {d^ - lA^f = -d^ + id^A^ + 2iA^d^ + A\ (5.122) 

so that 



Trln[-(9^-zA^)2 + ao] 



Trln(-a' + (To) + G(0) l df^xA^^x 



d'^x I d'^x'[id^A^{x) + 2iA^{x)d^]G{x - x') 



X [td'A^ix') + 2iA^ix')d']Gix'-x), 



(5.123) 



where 



G{x) 



d'^p 



{2n] 



-e^^-^G(p), 



(5.124) 



with 



G{p) 



1 



+ CTq 



(5.125) 



The first term on the RHS of (5.123) contributes to the saddle-point approx- 



imation. The second term can be written in momentum space as 
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where 



25. 



25, 



{2ny 



{p)A^{p)A,{-p), 



(5.126) 



{2tiY 
d'^k 



G{k) 
1 



d'^k 
(27r)^ 
d'^k 



{2k-p)^{2k~p),G{k-p)G{k) 
{2k-p)^{2k-p), 



{27rYk^ + ao J {27tY [{k - py + + (To) 



(5.127) 



Next we set = ctq and employ some tricks of dimensional regularization. 



First, by applying the operator Md/dM to both sides of Eq. (B.3) from 
Appendix [B| we obtain 



d'^k 



1 



2M^ 



d'^k 



1 



(5.128) 



(27r)'^P + M2 d-2j (27r)'^(fc2 + M2)2- 
Second, we can decompose the second integral appearing in the second line 



of Eq. (5.127) into its transversal and longitudinal parts, i.e. 



d^k [2k - pU2k - p). _ ^^^^^ / _p_^\ ^^^^^p,p. 



J {27iY [{k - py + M^]{k^ + M^) '^^^yt^'^ p2 J 'v^^ p2 ■ 

(5.129) 

By taking the trace of both sides of the above equation on one hand and 
contracting with p^p^y on the other hand, we can determine Dt{p) and Di{p): 
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and 



4M2 



d — 1 id 
+ (p2 + 4M2) 



+ M2)2 



[(A;-p)2 + M2](P + M2) 



(5.130) 



4M2 



(5.131) 



d-2 Jk (k^ + M^y 
where we have used the shorthand notation for the dimensional integrals 
of Chap. [2| Eq. (2.37). Thus, we obtain, 



-Dt{p) + 



4M2 



d-2 Jk {k^ + M^Y_ 



+ [A(P)-A(P)]^ 



S(p) - 



p 

PfiPv 
p2 



(5.132) 



where 



S(p) 



d 



-AM' 



+ (p2 + 4M2) [ 
Jk 



k (A;2 + M2)2 
1 



[(A;-p)2 + M2](P + M2) 



(5.133) 
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For large distances or, equivalently, in the infrared (i.e., p small), we can 
write ^ const p^, such that in real spacetime this fluctuation correction 
generates a Maxwell term. In order to obtain this result, we have to expand 
[(A; — p)^ + M^]^^ up to second order in p. Thus, we have 



1 



1 



{k-p^ + M^ k^ + M^ (k^ + M^) 



_^ 2p • k _^ Pf,Pu{4:kf,ku - V) _^ 



(F + M2)3 



(5.134) 



The above expression has to inserted in Eq. (5.133). In order to perform the 



integrals with help of the results from Appendix |Bj some simplifications are 
needed. First of all, rotational invariance implies 



k^ky 



5 



(5.135) 



lk{k^ + M^)^ d Jk{k^ + MY' 
Furthermore, the integral on the RHS of the above equation can be rewritten 
as 



Ik (F + M2)4 (F + M2)3 ^ Jk (P + M2)4 ■ 

It should also be noted that 

Now we are ready to express ^(p) in terms of the integrals evaluated in 
Appendix [B] to obtain 



123 



E(p) 



d-1 
P 



Hd) + ^ h{d) 1— /4(c?) 



d 



d 



(d- l)(47r)'^/2 

r 4 - - 

3ci V 2 



p'M^-^ f d 
3(47r)«'/2 V 2 



r 2 



2(4-d) / d 
2 7 d \ 2 



(5.138) 



where in the simphfications repeated use of the identity r{z + 1) = zr{z) 
was made. Therefore, after replacing back M — yfa^, we obtain at large 
distances the following low-energy contribution to the effective Lagrangian: 



'Maxwell 



1 / 1 

4 \e, 







(5.139) 



where F^,, = d^A^ - d^^A^ and 



1 

'''' " 3(4^ V 2 



(5.140) 



We have C2 = l/(127r) and C3 = l/(247r) for li = 2 and = 3, respectively. 
Note that even if a Maxwell term were not be present in the Lagrangian (i.e., 
eo — > 00), it would be generated by fluctuations in the paramagnetic phase. 

By using M = ^/oq as a mass scale, we can define a dimensionless gauge 
coupling constant / — M'^~^e^, where 



l + CdNelM'^-^' 
Therefore, the RG /3 function for / is 



(5.141) 
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M 



df 



iA-d)f + iA-d)c,Nf. 



(5.142) 



dM 



Note that by expanding around d = 4 we obtain exactly the one-loop /3 
function for the Abelian Higgs model in the minimal subtraction scheme [S] , 
although we have considered N large. 

We see that within the large approach the presence of a bare Maxwell 
term does not spoil the quantum critical point. Indeed, the behavior of 
for Cq large is the same as for M small (or g ^ Qc), provided 2 < < 4. 
However, as we have already mentioned, for the CP^ (i.e., for N = 2) model 
with a Maxwell term there is strong numerical evidence against the existence 
of a quantum critical point [HU EH]. However, there are indications that a 
doped version of the model exhibits quantum criticality [H]. 

Finite temperature at criticality 

The finite temperature analysis at large is similar to the one made for the 
non-linear a model, since for the leading order in 1/A^ we can set the gauge 
field to zero. This leads again to a thermal gap 



Let us next calculate the uniform susceptibility Xu at the QCP. In the case 
of the 0(3) non-linear a model it is obtained by considering the response to 
an external field twisting in the time direction of the direction field n. Thus, 
the following replacement holds 




(5.143) 
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{driif {dr-n -iHx nf; (5.144) 

see for example Ref. [6]. The CP^~^ model for = 2 is equivalent to the 
0(3) non-linear a model. In this case we replace 



drZa drZa - (5.145) 

for a field H = Hgz- We can generalize the above to the CP^^^ model by 
considering the generators of SU{N) instead of the Pauli matrices. Thus, at 
leading order the large free energy density at the presence of H becomes 
for D = 2 



n=— oo ^ ' 

+ \n{ujl + HuJn+p^ + Hy4 + m^)]- -^m^. (5.146) 



The uniform susceptibility is then defined by [HI EZj 

1 d'f 



(5.147) 



Therefore, 



2 



2^f-^^7 {27ry ool + + ^^^J (271)2(^2+^2 + ^2)2 

(5.148) 
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The Matsubara sums and the integrals are straightforwardly done to obtain 



m 



Air 



_ 1 



Tln(l - e 



-m/T\ 



(5.149) 



Upon substituting m{T) from Eq. (5.143), we obtain 



(5.150) 



5.3.4 Quantum electrodynamics in 2 + 1 dimensions 



We have seen in Section 15.11 that for f/ ^ t the Hamiltonian of the Hub- 
bard model becomes the Heisenberg Hamiltonian for an antiferromagnet 



[Eq. (5.10)]. In this Section we will consider a generalization of it in which 
the SU{2) symmetry is enlarged to SU{N). The aim of this generalization 
is to provide effective field-theoretic models that can be studied systemati- 
cally within a 1/A^ expansion. This approach was pioneered by Affleck and 
Marston [921 ES] long time ago with the physical motivation of better under- 
standing the Mott-insulating physics underlying the cuprates superconduc- 
tors. 

The SU{N) generalization of the Heisenberg model is given, up to irrel- 
evant constant terms by 



J 



(5.151) 



(id) 



where from now on summation over repeated Greek indices is implied, J = 
At^ /U [see Eq. (5.10)], and we have rescaled J — )■ J/N m. order to facili- 
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tate the large N approach. In addition, the fermions must satisfy the local 
constraint: 

flfia = N/2. (5.152) 

By performing a Hubbard-Stratonovich transformation, we can rewrite 
the Hamiltonian as 

H = jlXijf -Y.^vflfj<. + h.c. (5.153) 
The saddle-point solution at large N is given by the so called tt-Aux phase. 



Xi3 = Xoe'^'^ 

J2 Op = (5.154) 

r={ijki} 

corresponding to a flux amount of tt around a plaquette V. This gives the 
spectrum at large N: 



Ek — 4xo-\/cos^ kx + cos^ ky. (5.155) 

The above spectrum is gapless at the points (±7r/2, ±7r/2). 

The large N saddle-point approximation we just described corresponds 
to a Mott insulator without any broken symmetries. Thus, the solution 
is paramagnetic and preserves the symmetries of the square lattice. The 
obtained solution constitutes one of the simplest examples of spin liquid. 
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Let us now consider the fluctuations around the saddle-point. First, we 
would like to show that the elementary excitations near the nodes of the large 



N spectrum (5.155) are Dirac fermions. Here we will closely follow Refs. 
and [95]. Let us define 



XS — Xj+ej:+ey,j+eyi = Xj+ey,j- 



(5.156) 



Since the square lattice is bipartite, we can describe it as being composed by 
two interpenetrating sublattices, A and B, similarly to the mean-field theory 



for the Hubbard model at half-filling discussed in Sect. 5.1 Therefore, using 



Eq. (5.16), we can write 



H = -J2J2{'>^^4+ex,aCj^ + Xlc]+ey,aCj-+^-C' 



NL 



J 



|2 I I |2 I I |2 I I i2\ 
Xl\ + 1X2] + 1X3] + 1X4] ) ■ 



(5.157) 



Therefore, the Hamiltonian can be rewritten as 



Ixl (cos fci + i COS A;2) 

Ixl (cos /ci — i cos /i:2) 



(5.158) 
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where, as before in Sect. 5.1), the prime on the sum imphes that the sum is 
carried over the half of the Brilouin zone. Let us denote ciko- (ciko-) as the 
fermionic operator Cko- (cko-) near the node (7r/2,7r/2). Similarly, C2kcr (c2kcr) 
will denote the fermionic operator Cko- (cko-) near the node (— 7r/2, 7r/2). Thus, 
the Hamiltonian becomes approximately 



H 



k 



'^Ikcr '^Iko- '^2k(T '^2ko- 





(Tl 






-ai 



+ 1X1^2 



as 
cr2 



Ciko- 
Ciko- 
C2ko 
C2ko 



(5.159) 



where cti and a2 are the usual Pauli matrices. The above Hamiltonian is 
the one of four-component Dirac fermions in two spatial dimensions. If in 
addition we include the phase fluctuations of the link field Xij^ 



X/™ = Xe*^-, (5.160) 

the Dirac fermions will get gauged. The resulting effective theory that 
emerges in the continuum is quantum electrodynamics in 2 + 1 dimensions 
(QED2+1), with a Lagrangian given by 

N 

JO-f = Y,M^ + '('oA)^», (5.161) 

a=l 
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where we have introduced a bare "electric charge" cq to set the energy scale, 
and the standard notation (ji = 7^a^ is used. The Dirac matrices are in this 
case given by 



7o 



\ -(73 J 



72 



7i 



\ -(72 J 



V -a, J 

The gamma matrices above satisfy the anticommutation relation 



(5.162) 



7m7^ + lu'^i, = 26^^I, 



where J is a 4 x 4 identity matrix. The above relation implies 



(5.163) 



tr(7M7^^) = 45, 



(5.164) 



It is well known that massless QED2+1 featuring four-component Dirac 
spinors has a chiral symmetry [96j, in contrast with QED2+1 involving two- 
component spinors, which even admits a mass term breaking parity sym- 
metry. Indeed, for the QED2+1 above there are two 75-like matrices which 
anticommute with all the other Dirac matrices, namely. 



73 



/ 



/ 



75 



V 



/ 



-/ 



(5.165) 



with / being a 2 x 2 identity matrix. The chiral symmetry of the Lagrangian 
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(5.161 ) is then 



i) e^^^'^V- (5.166) 

Spontaneous breaking of the chiral symmetry leads to a dynamical mass 
generation in QED2+1 [98]. This mass is proportional to the chiral conden- 
sate, i.e., 

N 

m ~ (5.167) 

a=l 

The dynamical mass generation implies the spontaneous breaking of the 
SU{2) symmetry in the Heisenberg model, leading in this way to a stag- 
gered magnetization (for more details on this point, see Ref. [99j). A similar 
result |100j holds also in the case of effective QED theories for (i-wave super- 
conductivity pl)T|[Tn2]. 

The dynamical mass generation in QED2+1 is more easily obtained by 
studying the Schwinger-Dyson equations ^QSj . First of all, we need the photon 
propagator in the large N limit. This is obtained in a functional integral for- 
malism by performing the Gaussian integral over the Dirac fermions exactly. 
In this case Dirac fields are Grassmann variables (i.e., they anticommute), 
so that the effective action reads 

^cff = - A^Tr ln(^ + ieo4) . (5. 168) 

In order to obtain the photon propagator we simply expand the above effec- 
tive action up to quadratic order in A^. Although we are interested in the 
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d = 3 case, we will evaluate in d spacetime dimensions and set d = 3 

at the end. In momentum space we obtain in this way the quadratic form, 



where 



5'. 



eflF 



d'^p 
(27r)« 



{p)A^{p)A,{-p), 



with the fermionic propagator 



d'^k 



lfMGo{k)j^Go{p - k), 



(5.169) 



(5.170) 



Goip) = 

and a = Ne^. Gauge invariance implies p^T,^y{p) 
transverse, so we can write 



(5.171) 



0, so that T,^y{p) is 



T.,,{p)=p^n{p) (^5^ 

Taking the trace over the vector indices yields 



p2 



(5.172) 



{d - l)p^U{p) = a 



d'^k tr[7^|^7^(^ - . 
{27tY k'^{p-kf 



Next we use the identity 



(5.173) 



7mJ^7m = (2 - d)i), 



(5.174) 



and the trace formula (5.164) to rewrite Eq. (5.173) as 
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{d - ly n(p) = 4(c/ - 2)a 



Sk k ■ {p — k) 

{2Txf k^{p - ky 



Since in the numerator we have 



2 — k"^ — {p — ky 



k ■ p — k 



and assuming the rules of dimensional regularizationj^ we obtain 



(5.175) 



(5.176) 



(rf- l)n(p) = 2{d-2)a j 



d'^k 



(5.177) 



(27r)"'A;2(p- A;)2' 

The integral in the equation above was evaluated in Chap. [2| Sect. 2A_ It 



is the same as the one in Eq. (2.57), with the result of the integration given 



in Eq. (2.62). Therefore, 



n(p) 



2{d - 2)c{d) 

d-l ' 



(5.178) 



with c{d) given in Eq. (2.63). Thus, the gauge field propagator has the form 



DM 



p'^Y\.{p) 



p2 



(5.179) 



where Il{p) is given for = 3 by 



n(p) 



a 



l\p[ 



(5.180) 



Let us insert the propagator (5.179) in the one loop correction to the 
fermionic full propagator, 

^In writing the following equation we have used / cfikk^^ — J d'^k{p — fc)^^ = 0, which 
is true in dimensional regularization; see Appendix 151 
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G-'ip) = G,\p) + 



f2^ 



(5.181) 



with a dressed fermion propagator 



Gip) 



1 



— ifZij)) 

i^Z{p) + i:{p) ^ Z2(p)p2 + S2(p)' 



(5.182) 



Explicitly, we have, 



8 f d'k ^^[T.{k-p)+i{f-^)Z{k-p)]^^ 



N J {2Tif[Z^{k-p){k-pY + T.^{k-p)]\k 
8 [ d^k I/t[J:{k-p) + i{I/t-^)Z{k-p)]Iit 



N J (27r)3[Z2(fc-p)(A;-p)2 + S2(A;-p)]|A;|' 



(5.183) 



Thus, we obtain the selfconsistent equations: 



S(p) 



16 r d^k 



N J (27r)3 [Z^{k)k^ + S2(A;)] I A; + p\ 



(5.184) 



Z{p) = 1 - 



8 f d^k [e -p'^ + {k + pf]{k + p)-pZ{k) 



Np^ J (27r)3 [Z2(A;)A;2 + S2(A;)]|A; + p|3 



(5.185) 



In writing the above equations we have used the identity (5.174) and 



It^lt = 2p-kJ^- k^f. 



(5.186) 



After integrating over the angles, we obtain 
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S(p) 



dk 



dk 



kJ:{k){k + p- \k-p\) 

Z\k)k^ + T?{k) Z\k)k^ + S2(A;)_ 



(5.187) 



Z{p) = 1 



37r2A^ 



dk 



kZ{k) 



Z^{k)k'^ + Y?{k) p 



Z2(A;)A;2 + S2(fc)_ 

(5.188) 



From Eq. (5.187) we obtain 



p 



,dJ:{p) _ 8 
dp 



dk- 



eE{k) 



Ntt^ ./o Z^{k)k'^ + E'^{k)' 



(5.189) 



and therefore, 



d 
dp 



p 



dp 



8 P^S(p) 
'tc^N Z^{p)p^ + E^{p)' 



(5.190) 



A similar calculation using Eq. (5.188) yields 



d 
dp 



p 



, dZ{p) 
dp 



8 p'^Z{p) 
Z^{p)p^ + T.'^{p)' 



(5.191) 



In Refs. and [HH] Z{p) is considered as being approximately one, so that 



Eq. (5.191) is absent in their treatment. 



From Eq. (5.187) it follows the boundary condition 
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limpS(p) = 0. 



(5.192) 



Another important boundary condition is obtained from Eq. (5.187), since 



P- 



dp 



8 r ki:{k) 

^^^^ + iV7r2 '^^Z2(A;)P + S2(A;)' 



implying 



(5.193) 



P- 



dp 



p=a 



(5.194) 



which should be considered in addition to the boundary condition < S(0) < 



oo. Similarly, from Eq. (5.188) follows the boundary condition 



P- 



dZ{p) 



dp 



3[1-Z(a)], 



(5.195) 



while positivity of the spectral represention implies < Z{Q) < 1. 

In Ref. [96] it is found that a fermion mass is dynamically generated for 
all values of A^. Such a result can be obtained by considering the limit p — t- 



of Eq. (5.187) and replacing S(fc) and Z{k) in the integrand by their lowest 
order expansion in 1/A^, S(0) and unity, respectively. Doing this we obtain 
the gap equation 



E(0) 



8 S(0) 



''''P + S2(0)' 



(5.196) 



If we assume S(0) 7^ we can easily solve the above gap equation to obtain 



E(0) = ae 



(5.197) 
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reproducing the result of Ref. [96]. Under the same approximation used to 



derive Eq. (5.196), we obtain from Eq. (5.188) that 



Z{p) 



1 + 



8 



In 



v/p' + s2(o) 



a 



1 
3 



+ 



S2(0) S3(0) 



arctan 



P 



E(0) 



(5.198) 



If we insert in Eq. (5.198) the trivial solution S(0) = of the gap equation 



(5.196), we obtain 



Z{p) = 1 



37r2A^ 



In 



p\ 1 



a/ 3 



(5.199) 



The above equation implies the existence of an anomalous dimension rj de- 
fined through 



rj = — limp 



din Z{p) 



p^o dp 
and we obtain the known 1/A^ result |103] 



(5.200) 



On the other hand, the p — )• limit gives 



Z(0) = 1 



37r2A^ 



In 



S(0) 



a 



(5.201) 



(5.202) 



Substitution of the result (5.197) in (5.202) yields Z{0) = 2/3, which con- 
tradicts the value Z{0) = 3/4 obtained in Ref. [96] through an argument 
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involving the Ward identities The contradiction is removed if instead of 



Eq. (5.196) the following gap equation is used: 



S(0) 



8 S(0) 



dk 



k 



Z2(0)P + S2(0)' 



(5.203) 



In this case we obtain 



E(0) = aZ(0)e 



^Af7r2z2(0)/8 



(5.204) 



Z(0) = 1 + 



37r2ArZ(0) 



In 



aZ(0) 



implying that 



Z(0) 



(5.205) 



(5.206) 



Note that such a modification affects the value of the anomalous dimension, 
since in the limit of a vanishing gap we have 



Z{p) 



1 + 



37r2A^Z(0) 
32 



In 



97r2Ar 



In 



(5.207) 



instead of the standard 1/A^ result (5.199). Thus, we obtain instead of Eq 



(5.201) the result 



32 



97r2Ar 



(5.208) 



Thus, revisiting the earlier approach of Ref. [HS] taught us an impor- 
tant lesson: in order to not violate the Ward identities when linearizing the 



^This result is discussed in Ref. [19] of Ref. 
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problem, we have to account not only for S(0) but also for Z{0). 

We consider next an improved linearized problem given by the approxi- 
mated differential equations 



dp 
dp 



p 



p 



dp 

dZ{p) 
dp 



7r2A^Z2(0)p2 + S2(0)^ 



p'^Z{p) 



(5.209) 



(5.210) 



71^ Z2(0)p2 + S2(0)' 

The solutions obeying the initial conditions S'(0) = and Z'{0) = 0, while 
S(0) 7^ and Z{0) ^ 0, as required by the boundary conditions, are 



Z{p) = Z(0)2Fi 



S(p) = S(0)2Fi 



3 _ 3 3 3 5__Z2(0)p2 

4 " 4^' 4 ^4^' 2' S2(0) 



1 i 1 i 3 Z2(0)p2 
■7, 7 + 77: • 



4 4 "4 4 " 2' S2(0) J ' 
where 2F1 is a hypergeometric function and 



(5.211) 



(5.212) 



C 



32 



32 



97r2Z2(0)Ar' 



7 



7r2Z2(0)A^ 



1. 



(5.213) 



The solution (5.212) with Z{0) = 1 was obtained before in Ref. |104j . 

The CSB is more easily analysed in the regime where Z'^{0)p'^ ^ S^(0), 
in which case the above solutions simplify to 



Z{p) ^ ^-^Z{0)B 



mp 
m 



-3/2 



cosh — In 
2 



mp 
m. 



+ V>, (5.214) 



W,)~Mj![EMeosl^ln \^ 
^^P^^ 4 V Z(0)p'°n2 Ls(0)J 



(5.215) 
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where 



A 



± 



r(±C)(i±C) 



= 71^+^-1, c 



r(i7/2)(l + i7) 

(5.216) 



and 



= arccos 



C + C* 
2\C\ 



09 = - In 
^ 2 



A, 



(5.217) 



The boundary conditions can be used to determine S(0) and simphfy the 
above equations. For S(0) we obtain 



S(0) = aZ{0) exp < — arccos 



7 



nZ{0) IN ^ . ^ 

Y (sm^ + 7COS( 



(5.218) 



Note that the critical number of flavors is now modified due to the wave 
function renormalization at p = and given by demanding that 7 vanishes, 
i.e., 



N, 



32 



ch 



7r2Z2(0)' 



(5.219) 



Therefore, Eq. (5.218) can be rewritten as 



E(0) = q;Z'(0) exp <^ arccos 



N 



N, 



(sin 6' + 7 cos ( 



ch 



(5.220) 



Note that 6 is non-singular as N ^ Nch, approaching the value (2m + l)7r/2. 



with m = 1, 2, . . . . Thus, near = A^^h Eq- (5.220) becomes 



S(0) = aZ(0)exp 



27rn 

7 



(5.221) 



with 77, = 1, 2, ... , which agrees with Ref. [OE] for Z{0) = 1. Interestingly, 
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the functional form of the generated mass gap S(0) resembles the one in Eq. 



(3.63) for the KT transition. In the present case the number of fermionic 



degrees of freedom is playing the role of the temperature. 



Inserting (5.221) with n = 1 in (5.214) yields 



8 \a 



cosh 



2 \a/ 7 



. (5.222) 



Using Eq. (5.195) this can be further rewritten as 



Z{p) 



2{p/a) 



-3/2 



cosh(37r(^/7 + ^) + C sinh(37r^/7 + Lp) 



cosh 



2 \a/ 7 



(5.223) 

For > A^^ch the mass gap S(0) vanishes and the above equation becomes 
simply 



Z{p) 



p 

which defines the anomalous dimension 



(5.224) 



(5.225) 



Precisely at = A'^ch we have 1] = —0.08. On the other hand, for A^ ^ A'ch 
we obtain the following large A^ result: 



ch 



UN 3tt^Z^{0)N 



(5.226) 



The result (5.225 ) can be more easily obtained by considering Eq. (5.191 ) 
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at the critical point, in which case it becomes 



d 

dp 



P 



dp 



p" 



irmzipY 



(5.227) 



Eq. (5.227) can be rewritten as a set of two RG-hke equations: 



d\nZ{p) _ 
P IZ = -7(P), 



P- 



drf 



dp 



7(7-3) 



8 



(5.228) 



(5.229) 



dp TT^NZ^' 

Thus, rj = hmp_>o7(p)- This hmit is attained at the fixed point, which in 
this case is given by 



7± 



:i±c). 



(5.230) 



Only the solution 7_ makes sense since 7+ > 3, which would lead to a 
vanishing of the fermion propagator as p — t- 0. Therefore, we obtain that 



= 7_, in agreement with Eq. (5.225). 



Within the present calculation the determination of Z{0) is not so straight- 
forward, as the approximation we have used implies an anomalous dimension. 
On the other hand, since the anomalous dimension arises only for vanishing 
gap, it is reasonable to assume that rj coincides with the anomalous dimen- 



sion we have computed before, Eq. (5.208). Thus, we demand that rj = rj 



such that Zifi) = \/3/2 is obtained, contrasting with the Z{0) obtained be- 
fore, which is given by the square of the present result. Therefore, we obtain 



from Eq. (5.219) the critical value Nch = 128/(37r ) ^ 4.32, in full agreement 



with the result of Ref. |105j and the one of Ref. |106j in the case where the 
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anomalous dimension is computed up to order 

The physical interpretation of the above results is as follows. The num- 
ber Nch separates two different insulating regimes, namely, a magnetically 
ordered one at < Nch and another one having no broken symmetries 
at > Nch- The latter can be identified with a U{1) spin liquid state. 
The physically interesting case is = 2, which corresponds to an SU (2) 
Heisenberg antiferromagnet. Thus, the above Schwinger- Dyson calculation 
indicates that there is no spin liquid state in this case. However, there is 
numerical evidence |107] that the Schwinger-Dyson approach overestimates 
the value of N^ii- The numerical analysis made in Ref. |107] indicates that 
A^ch ~ 1-5. This result is in agreement with an earlier prediction |108j based 
on mathematical inequalities and symmetry arguments. If this result is cor- 
rect, it would imply that the two-dimensional Heisenberg antiferromagnet 
has a spin liquid ground state. 
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Appendix A 

The surface of a sphere in d 
dimensions 

In this appendix we will derive the formula for the surface of a sphere whose 
radius is unity in d dimensions. Such a factor occurs in many calculations in 
the text and we will derive it here as a simple exercise involving d-dimensional 
integrals. The trick is to consider the o?-dimensional Gauss integral 




(A.l) 



in two ways. In the first way we simply evaluate the above integral directly, 
since it is just an ordinary one-dimensional Gauss integral to the power d. 
Thus, 

/ = (A. 2) 
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Next, we use spherical coordinates in d dimensions and formally perform 
the integral over the d — 1 angular variables. We call the result of such an 
integration S^, which is precisely the surface of a unit sphere in d dimensions. 
In this way, we can also write, 

POO 

I = Sd drr^-'e-'\ (A. 3) 

Jo 

Now we make the change of variables u = r"^, such that. 







2 

SdTid/2) 



(A.4) 



where we have used the definition of the gamma function. From Eq. (A. 2) 
we obtain immediately, 



which is the desired result. 



149 



Appendix B 



Dimensional regularization and 
evaluation of some simple 
integrals in d dimensions 



In this Appendix we evaluate the following simple d-dimensional integrals 
and, at the same time, introduce the method of dimensional regularization 
of integrals: 



I d k 1 - , 



and 



' (B.2) 



(27r)<i(A;2 + M2)2- 

Both integrals can be obviously be done in 2 < d < 4 if an ultraviolet cutoff 
A is introduced. In the limit A ^ oo divergencies appear. Dimensional 
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regularization converts logarithmic divergencies due to the infinite cutoff into 
poles for some values of d. Moreover, divergencies involving positive powers 
of A "disappear". In order to understand the consistence of the technique, let 
us consider the integral Ji as a function of M and d, while demanding that it 
should be cutoff independent. This can be done by considering d somewhere 
in the complex plane, such that the integral is convergent. The result of 
the integration, as a function of can then be analytically continued in a 
larger domain (for more details, see Refs. [S] and ^). In this way, simple 
dimensional analysis yields 

/i(M,d) = ci(d)M"'-2, (B.3) 
where c{d) is a dimensionless function of d. We obtain similarly for /2, 

h{M,d) = C2{d)M''-\ (B.4) 

Note that if it is assumed that no cutoff is needed, the only scale available 
is M. This way of doing things is not always useful and sometimes we may 



need the explicit cutoff dependence. Note that Eq. (B.3) implies that Ji 
vanishes for all c? > 2 if M = 0, which is certainly not true if we regularize 
the integral with a cutoff. Indeed, if a cutoff is introduced, we have 

/i(0,rf) = ^3y, (B.5) 

thus leading to a positive power of A when d > 2. Thus, by demanding that 
/i(M, d) is independent of A for all d, we obtain that this condition actually 
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imposes 



d^k 1 ^ , 
= 0. (B.6) 



The above equation is one of the properties required by dimensional regular 
ization. In fact, an even more general formula should be used, namely. 



/ 



' 0, (B.T) 



where a e C. 

We now proceed with the explicit evaluation of the constants Ci(d) and 
C2{d). To this end, note that the integral h can be written as 

The d-dimensional Gauss integral is easily performed to obtain. 

By changing the variable as m = M^A, we obtain, 

h = 7^ r duu-'/'e-^ = (l-i), (B.IO) 

where we have used the integral definition of the gamma function. 

The integral I2 is obtained from Ji via simple differentiation of Ii with 
respect to M^. We obtain. 
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Note that the integral Ji has poles for = 2 and for = 4, while I2 has 
only a pole at c? = 4. 

Now we will consider the most general case, i.e., we will evaluate the 
integral 



d'^k 



(B.12) 



Using the identity (2.43), we have 



r(a)7o J i^^r 



1 



00 



(47r)d/2r(a 

l\/fd—2a foo 

ikf^-^" f d\ , , 

r a-- , (B.13) 



(4vr)<i/2r(a) V 2 

where from the first to the second line we have performed the d-dimensional 
Gauss integral in k and from the second to the third line we made the sub- 
stitution u = M^A. From the third to the last line we used the definition of 
the gamma function. 



Of course, for a = 1 and a = 2 the integral (B.13) reproduces the results 
Ii and I2, respectively. 
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Appendix C 

Evaluation of the integral over 
a Bose-Einstein distribution in 
d dimensions 

Let us consider the following integral: 

^^^^"'^^ - J {2nr e^\^\^ - 1 

where c > and 5*^ is given in Appendix |Aj The exponent z > defines the 
power law for the spectrum of the corresponding Bose particle. For example, 
a free particle will have z = 2, while phonons have z = 1. 

After performing the change of variables u = ck^, the integral becomes, 
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In order to calculate the integral appearing in the expression above, we con- 
sider the integral representation of the gamma function, 



oo 

s-1 -t 



T{s) = / dtf-'e-\ (C.3) 







and make the replacement t ^ nt, where n is a positive integer. This leads 
to 



- ' dtt^-^e-"*. (C.4) 



n" 



poo 



Now we sum over n from 1 to oo, 

CXD ^ n OO 

n=l n=l 

On the LHS of the above equation appears the definition of the zeta function, 

oo ^ 

n=l 

while on the RHS features a trivial geometric sum. 



1 1 

z-^ ^ — p-t pt — ^ ^ 

Therefore, 



1 - e-* e* - 1 

n=l 



oo ^.-1 

e^- 1 



C(s)r(.) = / di-^. (C.8) 
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Using the above result in Eq. (C.2), we obtain finally. 



156 



Appendix D 



Matsubara sums 



The computation of Matsubara sums is a standard method in both relativistic 
and non-relativistic quantum field theory; see for example Ref. |109] . 

Let us evaluate here two bosonic Matsubara sums frequently used in the 
text, namaly, 



where Un = 27m/ (3. 

Although, strictly speaking, the first sum does not converge, it can nev- 
ertheless be done using regularization methods. Moreover, as we will see, it 
is related to the second sum, which is convergent. 

Since 5*2 is convergent, let us calculate this one first. 




(D.l) 



and 




(D.2) 
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Figure D.l: Initial contour integration in complex plane used for performing 
the Matsubara sum 5*2 . The red crosses indicate the position of the Matsub- 
ara modes. 



In order to calculate 5*2, we rewrite the sum in terms of a contour integral 
over an appropriate contour in such a way as to set the sum of the residues 
at poles given by the Matsubara frequency modes. The two straight lines 



shown in Fig. D.l enclose all Matsubara modes. We can imagine a finite 
number of modes with a closed curve around all of them. In the limit of 
an infinite number of modes, the closed curve becomes two parallel lines 
having the Matsubara modes between them. This picture explains why the 
arrows are pointing the way they do: if to the closed curve a counterclockwise 
orientation is given, in the limit of an infinite number of Matsubara modes 
the resulting straight line on the right will point upwards, while the one on 
the left will point downwards. Thus, the sum S2 is equivalent to 
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)( 

X 
X 



Figure D.2: New contour integration in complex plane used for performing 
the Matsubara sum 5*2 . The contours C'l and C'2 enclose the poles z = ±e. 



Jc,uC2 27r« - e2 e^^ - 1 
The contour Ci U C2 encloses the poles of the function l/{e^^ ~ 1); but not 
the two poles of — e^). The straight lines Ci and C2 can be deformed 



into the curves C[ and C2 shown in Fig. D.2 without encountering any 
singularity. Note that the deformed contours go clockwise around the poles 
z = ±e, which introduces an extra minus sign. Thus, the residue theorem 
yields, 



82 = ^ coth ] . (D.4) 



2e V 2 

Now we will consider the sum Si. In contrast to the case of 5*2, this sum 
is related to only one pole at ^ = e in the contour integration. The initial 



contour is shown in Fig. D.3 Now by deforming the contours Ci and C2 
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Figure D.3: Initial contour integration in the complex plane used for per- 
forming the Matsubara sum Si. 



into the new contours C[ and illustrated in Fig. D.4, we see that does 



not enclose any pole and thus contributes nothing in the evaluation of the 



sum. Therefore, we obtain, 



^1 = (D.5) 



The sum Si can be related to S2 in the following way. 



^1 



^ 00 

E 



e + iu„ 



/3 



n ' 



1 



n ' 



(D.6) 



The second term in the second line of the above equation can be thought to 
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vanish, as to every positive term there is a corresponding negative term of 
identical magnitude. Thus, it seems that we just have to plug the result of 
the sum 5*2 into the above result and we are done. However, this is not quite 
correct. The point is that S2 involves two poles, z = ±e, corresponding to a 
"particle" and an "anti-particle" mode, while 5*1 involves just a single pole 



a.t z = e. Actually the second sum in Eq. (D.6) is not convergent. This is 
consistent with the the already mentioned fact, that the sum (D.l) is not 
convergent. Since, 



Si = - coth ( ^ ) + ^3 = - + — + ^3, (D.7) 

where 



1 \ ^ 100. 



n=—oo 



and using Eq. (D.5), we obtain that 



^3 = -l- (D.9) 

The factor 1/2 subtracted via the sum 5*3 reflects the "vacuum" con- 
tribution associated to the particle-antiparticle pair. In 5*2 this leads to a 
finite zero-temperature contribution, which is absent in the non-relativistic 
Matsubara sum. 
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Figure D.4: Contour integration in complex plane used for performing the 
Matsubara sum 5*1. The contour C[ encloses the pole z = e, while no pole is 
enclosed by the contour C!^- 
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Appendix E 



Classical limit of the 



polarization bubble of the 
dilute Bose gas 



After setting Xq = fi Eq. (4.44) can be rewritten as 



n(za;,p) = 4m / 77— t;tT1b ' 



(2tiY \2m J \2miu — p"^ — 2p ■ q 2mza; + — 2p ■ q 

(E.l) 

In the classical approximation we write nb{x) ~ T/x and the polarization 
bubble can be rewritten as 

U{iuj, p) = 4m^T{I+ - /_), (E.2) 

where 
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^ J (27r)'^2ma; + i(2p-q±p2)g2- ^ ' > 

The integrals I± can be evaluated using the Feynman parameters [8], 



± 



/•CO i-oo 

Jo Jo (2vr)'^ 



After evaluating the Gaussian integral over q we obtain 



d/2 



W Jo 

T^§^r(rf/2 - i)r(3 - rf)/- ( 1 T 



g Xi{2mLU±ip^)^-p^\l/X2 



(E.5) 



Substituting the above expression back into (E.2) we obtain Eq. (4.61). 
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